12

2.3. Представление геологической модели

В системе DV-Geo представление модели геологического объекта имеет три стороны – физическую реализацию (каталог файлов), интерфейс в программной среде (окно-список объектов) и связующий файл (динамический документ). Названные стороны образуют взаимосвязанное и согласованное единство, каждая из них безусловно необходима, и ничего сверх названного в представление модели не входит. Соответственно, все без исключения элементы модели (скважины, поверхности, каротажные кривые и т.д.) в обязательном порядке присутствуют в составе каждой из трех названных сторон - как файлы в определенных каталогах, как имена и пиктограммы в окне-списке объектов и как строки-описания в связующем файле. Рассмотрим три названные стороны представления модели более подробно.

2.3.1. Каталог файлов

Первая сторона представления – физическая реализация модели. Это структура каталогов с файлами, имеющими отношение к модели конкретного геологического объекта (фиг. 4).

Данная структура каталогов создается автоматически там и тогда, где и когда пользователь сохраняет модель (или проект) под новым именем. При этом создаваемый каталог верхнего уровня получает имя проекта, вложенный в него каталог получает фиксированное имя DEFAULT, и в нем, в свою очередь, создаются вложенные каталоги с фиксированными именами так, как показано на фиг. 4. В данном случае в рабочем каталоге «D:\DvGeo\GeoDat» был создан проект под именем «cge». В начальный момент почти все вложенные каталоги пусты. Их наполнение производится системой по мере включения в модель новых объектов, что происходит в результате или загрузки внешних данных, или выполнения расчетных процедур, или создания элементов модели непосредственно пользователем. Появляющиеся при этом файлы являются принятой в системе реализацией скважин, каротажных кривых, стратиграфических отметок, поверхностей, карт, кубов параметров и многих других объектов, в совокупности составляющих геологическую модель. Хотя значительная часть этих файлов имеет текстовый формат и доступное пониманию содержание, пользователь не должен работать с ними впрямую. Система DV-Geo предоставляет доступ к этим файлам через специальный интерфейс – окно «Список объектов».

13

Фиг. 4. Структура каталогов с файлами реализации геологической модели

2.3.2. Окно «Список объектов»

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

Окно «Список объектов» принадлежит к особому оконному классу, специально разработанному для отображения древовидных списков. Объекты в таком окне распределены по группам, образующим иерархическую структуру. Наименования групп верхнего уровня

14

фиксированы, каждая такая группа предназначена для размещения объектов определенных типов. Например, группа «SKV» служит для размещения скважин и данных ГИС, «Correlation» - для скважинных отметок, «Region» - для регионов, поверхностей и карт, «Pack» - для кубов параметров и т.д. Исключением является группа «Tools» - ее объектами являются процедуры, и о ней будет рассказано особо.

Каждая группа в окне «Список объектов» может содержать вложенные группы, созданные пользователем для упорядочивания объектов. Находящиеся в группах объекты исходных типов очень часто одновременно являются группами для своих дочерних объектов. Например, каждый регион является группой для заданных на нем поверхностей, поверхность является группой для производных карт, карта является группой для наборов изолиний и т.д. Чтобы раскрыть (или, наоборот, скрыть) содержание группы необходимо дважды щелкнуть на ней левой кнопкой «мыши».

Окно «Список объектов» при работе с моделью всегда только одно, хотя окон отображения модели (3D, 2D и других) может быть сколько угодно.

Фиг. 5. Окно «Список объектов»

2.3.3. Связующий dvg-файл

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

15

считывается, когда пользователь выбирает в главном меню системы DV-Geo команду «Файл\Загрузить» и сохраняется по команде «Файл\Сохранить» или «Файл\Сохранить как» (одновременно с сохранением или обновлением содержания описанной выше структуры каталогов «GeoDat\Имя проекта\ …»). Dvg-файл имеет текстовый вид и его можно просматривать с помощью программы Notepad или WordPad. Большую его часть составляют описания в установленном формате элементов, отображаемых в окне-списке – самих объектов и групп, по которым они распределены. В описании каждого объекта присутствует его имя, автоматически присвоенный ему уникальный номер, ссылка на файл, в котором этот объект реализован, а также ссылка на группу, в которую он включен.

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

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

Источник описания представления модели - «Руководство пользователя», дополнительный том 1 «Организация данных и процедур», раздел 1.1.

2.4. Организация процедурной части DV-Geo

Исходный код системы DV-Geo состоит из двух частей. Одна часть исходного кода, выполненная на языке С++, от пользователя закрыта. Эта часть включает ядро системы, в которое входят, например, процедуры визуализации объектов модели. Закрытость этого кода гарантирует, в частности, надежную работу упомянутой визуализации. Пользователю, кстати говоря, от этой части системы ничего, кроме надежной работы, и не нужно.

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

В DV-Geo названная проблема решается следующим образом:

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

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

16

процедур на языке Th в состав системы (в упомянутую выше группу «Tools» окна «Список объектов», а также в меню объектов разных типов, представленных в том же окне);

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

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

Организация процедурной части системы DV-Geo почти точно повторяет описанную выше организацию модели, т.е. также имеет три стороны. Первой стороной (физической реализацией) является каталог THLIB с текстами процедур. Второй стороной (интерфейсом в программной среде) является процедурная группа «Tools» окна «Список объектов» и многочисленные меню самих объектов. Третьей стороной является связующий файл, в котором названия процедур в окне-списке и в меню объектов ассоциируются с текстами процедур. Связующим является файл «DvGeo\thlib\rout\thtolls.tht».

Пользователь может включать в группу «Tools» и в меню объектов новые процедуры, а также назначать для существующих названий процедур и пунктов меню другие файлы с текстами процедур. В ходе диалога по поводу производимых изменений ему необходимо указать, должны ли эти изменения вноситься в библиотеку (опция LIB), или они должны касаться только загруженной модели (опция USER). В первом случае модифицируется связующий файл «thtolls.tht», т.е. фактически изменяется состав системы DV-Geo на данной платформе. Однако сохранение этих изменений за пределами текущего рабочего сеанса происходит лишь тогда, когда при запуске системы указывается специальный ключ. Во втором случае изменения фиксируются в файле «…\Имя модели\anything.bin» и действуют при работе только с данной моделью. Право на такие изменения имеет каждый пользователь.

Источник описания организации процедур - «Руководство пользователя», дополнительный том 1 «Организация данных и процедур», раздел 1.2.

После того, как мы познакомились с принятой в DV-Geo организацией данных и процедур, мы можем приступить к практическому освоению системы. Это освоение мы проведем в два этапа. На первом этапе мы будем изучать готовую модель. Это реальный проект, представленный в инсталляционном наборе DV-Geo под условным названием «CGE». В ходе этой работы мы освоим пользовательский интерфейс системы, а также увидим, что собой представляет построенная геологическая модель. На втором этапе мы пройдем весь путь построения той же самой модели от начала и до конца.

17

Часть I. Знакомство с геологической моделью «CGE»

Из раздела «Введение» мы знаем, что геологическая модель есть не что иное, как набор файлов. Поскольку здесь уже не «Введение», скажем точнее. Геологическая модель есть некоторый набор исходных и производных абстракций (траекторий скважин, каротажных кривых, корреляционных отметок, интерполированных поверхностей, кубов параметров и других), каждая из которых представлена форматированным набором данных. Суть моделирования состоит в том, что образы этих абстракций и вычисления на их основе рождают новые и вполне конкретные знания - о генезисе геологического объекта, о свойствах резервуара, об объеме запасов, о возможностях их извлечения и т.д.

Наше рассмотрение готовой модели будет состоять в том, что сначала мы познакомимся с полным списком элементов модели в окне «Список объектов» (глава 1) , а затем с отображением этих элементов в окнах визуализации четырех основных типов (главы 2-5). Для лучшего понимания, что это такое, используем следующую аналогию. Список объектов – это модельная реальность, она одна. Окна визуализации - это телевизионные камеры нескольких типов, отображающие модельную реальность. Этих «камер» может быть сколько угодно. Пользователь сам расставляет их там, где считает нужным. За пределы аналогии выходит лишь то обстоятельство, что через окна визуализации мы можем не только наблюдать, но и воздействовать на модельную реальность.

18

1. Представление модели в окне «Список объектов»

1.1. Главное окно системы DV-Geo. Открытие проекта

Для запуска системы дважды щелкните левой клавишей мыши по соответствующему ярлыку на рабочем столе. После этого откроется главное окно DV-Geo, которое имеет следующий вид

(фиг. 6):

Фиг. 6. Главное окно системы DV-Geo

Имя «Undefined» в заголовке главного окна указывает на то, что пока еще никакая модель не загружена. Из этого состояния мы можем начать создание новой модели, имя которой мы должны будем задать при первом сохранении. В правом верхнем углу главного окна

расположены кнопки , предназначенные для сворачивания окна, его максимизации и выхода из программы. Эти же функции доступны через системное меню главного окна,

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

Загрузим в DV-Geo проект «CGE». Для этого активируем команду главного меню «Файл\Загрузить» и выберем расположенный в рабочем каталоге GeoDat файл «cge.dvg». В строке «Тип файла» должно быть предварительно указано «Проекты DV-Geo» (фиг. 7).

19

Фиг. 7. Выбор в рабочем каталоге файла «cge.dvg»

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

Если dvg-файл случайно испорчен, его чтение прерывается и выводится сообщение с указанием номера проблемной строки. В этой ситуации необходимо просмотреть dvg-файл в текстовом редакторе и найти соответствующую строку. Как правило, проблему можно понять и исправить. В виде исключения, осторожное редактирование dvg-файла допускается. Второй способ выхода из названной ситуации – переименовать и заново использовать bak-версию dvg-файла (его вид, предшествующий последнему изменению), которая сохраняется в папке с именем проекта при каждом выполнении команды «Сохранить проект».

После успешного окончания загрузки имя dvg-файла (включая путь) появляется в заголовке главного окна. В строке сообщений отображается текст «Модель “…\DvGeo\GeoDat\cge.dvg” загружена». Одновременно открывается окно «Список объектов» (фиг. 8).

20

Фиг. 8. Результат успешной загрузки модели

Скорее всего, начальное положение и размеры окна «Список объектов» у Вас будут иными. Причина этого в том, что DV-Geo запоминает последние экранные координаты всех своих окон, и затем воспроизводит их в ходе очередного сеанса. Если пользователь имел накануне два монитора, а сейчас имеет только один (или если он перед этим использовал видеорежим с более высоким разрешением), то некоторые окна он может не увидеть, поскольку их координаты будут за границами экрана. В этом случае необходимо удалить из каталога GeoDat файл winpos.txt. После удаления названный файл будет создан заново, и все окна будут открываться в левом верхнем углу экрана с минимальными начальными размерами.

1.2. Окно «Список объектов». Группы объектов

Окно «Список объектов» является главным окном проекта. О его назначении мы уже говорили во Введении, в разделе 2.3. Оно представляет собой интерфейс базы данных проекта. Тип этого окна – древовидный список. Окно «Список объектов» в проекте всегда только одно.

Параллельно открыть второй проект (второе полнофункциональное окно «Список объектов») в DV-Geo нельзя. Однако можно параллельно открыть так называемый «Общий список» другого проекта. Это окно такого же типа, в которое может быть занесен ряд объектов другого проекта. Использование общих списков делает возможным обмен данными между проектами.

Хотя, как мы сказали, окно «Список объектов» только одно, можно параллельно открывать древовидные окна-списки с отображением отдельных ветвей иерархии объектов. Например,

21

можно открыть окно-список с отображением объектов одной скважины, или одного региона, или одной группы процедур, и.т.д. Такая возможность помогает быстрее находить нужные объекты. В отношении своих объектов окна-ветви функционально равноценны главному окну «Список объектов».

На фиг. 9 окно «Список объектов» показано вместе с пояснениями, касающимися содержания основных групп верхнего уровня. Повторим, открывать и закрывать группы можно двойным щелчком левой кнопки мыши. Щелчок правой кнопкой мыши позволяет открыть меню, ассоциированное с соответствующим объектом или группой.

Фиг. 9. Список объектов с пояснениями

Исключая группу процедур, окно «Список объектов» показывает нам всю совокупность элементов одной геологической модели. Ничего сверх этого списка в модели нет. Сейчас мы рассмотрим основные группы данного списка и узнаем, из каких смысловых абстракций (и выражающих их форматированных данных, файлов) состоит модель «CGE». Увидеть образы всех этих элементов мы сможем позже, когда будем знакомиться со средствами визуализации модели в окнах 3D и 2D.

1.3. Группа «Tools» (процедуры)

Как уже было сказано во Введении (в разделе 2.4), группа «Tools» является особой - она содержит процедуры на языке Th. В отличие от других групп, ее содержание является одинаковым для всех моделей (если только пользователь не увлекается написанием собственных процедур). Процедуры упорядочены по подгруппам. О назначении процедур

22

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

Фиг. 10. Меню, ассоциированное в «Списке объектов» с каждой процедурой

Пункт «Выполнить» приводит к запуску процедуры. Пункт «Редактировать» открывает окно редактора Th-процедур с загруженным текстом процедуры. Другая возможность открыть

редактор Th-процедур - нажать кнопку на панели главного окна. Пункт «Установить» открывает стандартный диалог, который показывает имя и размещение файла с текстом процедуры, ассоциированного с данным пунктом в списке «Tools». В ходе этого диалога пользователь может выбрать для данного пункта другой исходный файл. Использование редактора Th-процедур описано (на примере процедуры загрузки скважин) в дополнительном томе 1 «Организация данных и процедур», в разделе 3.1. Описание языка Th дано в дополнительном томе 2 «TH Процессор».

Несмотря на свою объемистость, группа Tools содержит лишь малую часть процедур DV-Geo. Большую часть составляют процедуры, ассоциированные с объектами других групп. Доступ к этим процедурам осуществляется через меню объектов. Как это сделано, описывается в следующем разделе.

1.4. Группа «SKV» (скважины, данные по скважинам)

Внутри группы «SKV» находятся две группы скважин - «Well» и «Well_1». Названия этих вложенных групп определяет пользователь при загрузке массивов скважин из внешних файлов (см. дополнительный том 1 «Организация данных и процедур, раздел 3.2»). Модель «CGE» построена только по скважинам из группы «Well», скважины из группы «Well_1» не используются.

Раскроем группу скважин «Well». Внутри нее находится множество групп, каждая из которых имеет имя скважины (фиг. 11). Один элемент в группе с именем скважины также имеет имя скважины – это ее траектория. Еще в группе с именем скважины обычно присутствуют группы «Gis» (совокупность каротажных кривых) и «Rgi»(совокупность интервальных разбивок). Имена этих групп определяются пользователем при загрузке каротажных кривых и интервальных разбивок. Групп данных ГИС и РГИ может быть несколько. Еще в группе с

23

именем скважины присутствует объект «Конструкция» (форма для описания конструкции скважины) и группа «Perforation» (данные о перфорации).

Фиг. 11. Состав данных и процедур группы «Well»

С группой «Well» ассоциировано меню, показанное на фиг. 11. Легко убедиться, что корневой пункт «Средства» присутствует в меню почти каждой группы или объекта. Его особенность в том, что за ним следуют процедуры, реализованные на языке Th. В меню пункта «Средства» всегда присутствует выделенный пункт «Определить средства». Он открывает окнодревовидный список «Th List», в котором группа верхнего уровня является аналогом рассмотренной выше группы «Tools». Отличие процедур меню от процедур группы «Tools» состоит в том, что первые всегда получают адрес объекта или группы, с которыми они ассоциированы. Изменяя состав процедур в окне «Th List» можно изменять состав меню пункта «Средства», то есть изменять меню объектов и групп.

Вернемся к группе с именем скважины. Меню, ассоциированное с этой группой, показано на фиг. 12. На первой панели меню рядом с пунктом «Средства» присутствуют пункты «Добавить керн» и «Редактировать\Инклинометрия». Пункты меню, не входящие в группу «Средства», реализованы на С++ и изменить их нельзя. Выбор пункта «Инклинометрия» открывает окно для просмотра и ручного редактирования таблицы инклинометрии (справа на фиг. 12). В проекте «CGE» траектории всех (загруженных на текущий момент) скважин представлены как вертикали, поэтому таблицы инклинометрии состоят только из двух строк. Для редактирования записей таблицы необходимо выделить соответствующую ячейку таблицы мышью и нажать «Enter» Для добавления строки после строки с выделенной ячейкой необходимо нажать клавиши «Ctrl» + «Ins». Удаление строки с выделенной ячейкой производится нажатием клавиш «Ctrl» + «Del».

24

Фиг. 12. Меню группы с именем скважины и таблица инклинометрии

На фиг. 13 показано содержание групп «Gis» и «Rgi», а также меню, ассоциированное с каждой каротажной кривой. На первой панели меню рядом с пунктом «Средства» присутствуют пункты «Редактировать» и «Информация». (как мы уже знаем, эти процедуры реализованы на С++). По первому пункту открывается окно для просмотра и ручного редактирования значений каротажной кривой (справа вверху), по второму – окно с общей информацией о каротажной кривой (справа внизу). Редактирование таблицы с данными каротажа производится аналогично редактированию таблицы инклинометрии.

Фиг. 13. Состав групп «Rgi» и «Gis», меню каротажной кривой «PS», таблица и окно «Информация» кривой «PS»

В каких каталогах система DV-Geo хранит загруженные в нее координаты скважин, данные инклинометрии, каротажные кривые и разбивки Rgi проекта «CGE»? В каких форматах она хранит эти данные?

25

Координаты (устьев или пластопересечений) скважин группы Well хранятся в текстовом файле CGE\default\LST\Well\xyz.txt.

Данные инклинометрии скважин группы Well хранятся в бинарном файле

CGE\default\LST\Well\wells.vfs.

Данные каротажа скважин группы Well хранятся в каталоге CGE\default\LAS\Well в бинарных файлах с расширением «lfs».

Интервальные разбивки скважин группы Well хранятся в в каталоге CGE\default\RGI\Well в бинарных файлах с расширением «rfs».

1.5. Группа «Correlation» (корреляционные маркеры)

Объектами группы «Correlation» являются корреляционные маркеры на скважинах. Один корреляционный маркер (или репер) представляет собой множество отметок одной и той же стратиграфической поверхности на всех (где она прослежена) скважинах.

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

Фиг.14. Корреляционные маркеры модели «CGE». Показано меню группы «Correlation» и меню одного маркера

Модель «CGE» содержит в группе «Correlation» семь маркеров (фиг. 14). Все они были созданы вручную в ходе корреляции данных ГИС в окне «Профиль» (см. ниже). Из названий маркеров следует, что они определяют кровлю и подошву стратиграфических горизонтов 71с

26

и 72с, водонефтяной контакт (OWC7_1), кровлю и подошву коллектора горизонта 7_1. Меню маркера также показано на фиг. 14. Пункт «Вычислить» открывает диалоговое окно процедуры автоматической корреляции, работа с которой описана в Томе 3-2. Пункт «Копировать» позволяет сделать копии маркеров с заданным смещением по кабельной длине или по абсолютной глубине. Пункты «Редактировать» и «Информация» открывают окна, показанные на фиг. 15. Редактирование таблицы маркера производится аналогично редактированию таблицы инклинометрии.

Фиг. 15. Окно для редактирования корреляционного маркера и окно «Информация»

Имя каталога хранения маркеров …\CGE\default\CORRL. Файлы маркеров имеют расширение «crr». Формат файлов текстовый.

Все остальные пункты меню маркера и группы «Correlation» реализованы на языке Th, то есть пользователь может самостоятельно разобрать их исходный код по строкам. Доступ к текстам процедур осуществляется через пункты меню «Определить средства». Через эти же пункты пользователь может наращивать показанные меню.

1.6. Группа «Region» (регионы, поверхности, карты)

Основное содержание группы «Region» - поверхности и карты. В группе присутствуют разнотипные объекты, имеющие свою иерархию (фиг.16). Объектом верхнего уровня является регион – прямоугольная область на плане X,Y, разбитая на ячейки регулярной сеткой. В модели «CGE» используется только один регион «1». Регионы создаются через меню группы

«Region».

27

Фиг. 16. Иерархия объектов в группе «Region». Показано меню объекта «Поверхность»

На регионах определяются поверхности и карты. Каждая поверхность создается через меню созданного ранее региона, что позволяет ей использовать его сетку. При этом она становится его дочерним объектом. Можно также создать поверхность на регионе через меню его поверхности, созданной ранее. В обоих случаях используется пункт «Добавить поверхность»

(фиг. 16).

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

Имея поверхность, можно создать уже ее дочерний объект - карту. Главное отличие между ними в том, что значения поверхности ассоциированы с узлами сетки на плане X, Y, а значения карты ассоциированы с центрами ячеек сетки на плане X, Y. Поэтому на одной и той же сетке карта имеет на единицу меньше значений вдоль X и вдоль Y, чем поверхность. Далее, через меню карты можно создать дочерний объект карты – «набор изолиний».

Регионы хранятся в каталоге …\CGE\default\CORRL. Файлы с описанием регионов имеют расширение «trg». Формат файлов текстовый.

Поверхности хранятся в каталоге …\CGE\default\CORRL\SURF. Файлы поверхностей имеют расширение «crs». Формат файлов текстовый.

Карты хранятся в каталоге …\CGE\default\CORRL\MAP. Файлы карт имеют расширение «crm». Формат файлов бинарный.

28

Изолинии хранятся в каталоге …\CGE\default\LEVLINE. Файлы изолиний имеют расширение «flv». Формат файлов бинарный.

1.7. Группа «Pack» (кубы параметров)

Объектами типа «Pack» являются кубы параметров. Для задания куба необходим регион (сетка на плане X,Y) и минимум две поверхности, заданные на этом регионе. Поверхности нужны для определения верхней и нижней границы «куба» (термин «куб» применяется здесь со значительной долей условности). Чаще при задании куба используют не две, а некоторую стопку предварительно подготовленных структурных (стратиграфических) поверхностей.

Для создания куба необходимо в меню группы «Pack» выбрать пункт «Добавить куб». Такой же пункт есть в меню каждого существующего куба. Во втором случае геометрия нового куба (регион, структурные поверхности, разбиение объема на ячейки) воспроизводит геометрию существующего. Пункт «Добавить куб» создает шаблон куба. Расчет значений в ячейках куба производится с помощью диалогового окна «Вычислить». Основной способ расчета кубов – интерполяция данных Gis или Rgi в скважинах. Состав группы «Pack» проекта «CGE» и меню, ассоциированные с кубами, показаны на фиг.17.

Фиг. 17. Состав группы «Pack» проекта «CGE» и меню, ассоциированные с кубами

Кубы хранятся в каталоге …\CGE\default\CORRL\PACK. Файлы кубов имеют расширение «crp». Формат файлов бинарный.

29

1.8. Группа «Area» (области)

Объектами типа «Area» являются замкнутые и незамкнутые контуры (линии) на плане X,Y. Они используются для отображения на плане границ лицензионных участков, контуров ВНК и ГНК, профилей, разрезов, линий разломов, областей с разными категориями запасов, рек, озер и т.д. Замкнутые контуры могут быть границами области расчета поверхности, куба, селекторами скважин. Незамкнутые линии служат образующими при построении разломов в пространстве 3D. При необходимости, в качестве объекта «Area» можно создать и сохранить любой график в координатах X,Y.

Состав группы «Area» модели «CGE» и меню, ассоциированные с группами объектов этого типа, показан на фиг. 18. В группе «Raz» находится линия разлома. Группа «kat» содержит контуры областей, в которых подсчет запасов производится по категориям B и С1. В группе «OWC» находятся контуры внешних и внутренних водонефтяных контактов. Кривые «L0» и «OWC_vn» использовались при построении контуров ВНК.

Фиг. 18. Группа «Area» модели «CGE» и ассоциированные меню

Контуры типа «Area» хранятся в каталоге …\CGE\default\AREA. Соответствующие файлы имеют расширение «aaa». Формат файлов текстовый.

1.9. Группа «Objects» (таблицы)

Объекты данного типа представляют собой таблицы вида «скважина – числовое значение (или дата)», а также «скважина – два числовых значения (или две даты)». Данный тип объектов разработан специально для разного рода фильтров, производящих отбор скважин по значению (или по двум значениям) некоторого параметра. Например, с помощью таблицы «Скважина – дата бурения» специализированный фильтр может визуализировать историю разбуривания площади. Названный тип данных очень близок к описанному ранее типу «корреляционный маркер», который также может использоваться для разделения скважин на подмножества.

30

Фиг. 19. Группа «Objects» модели «CGE» и ассоциированные меню

Таблицы типа «Objects» хранятся в каталоге …\CGE\default\PROPS. Соответствующие файлы имеют расширение «wdt». Формат файлов текстовый.

1.10. Группа «Lattice» (решетки)

Для гидродинамических симуляторов необходимо пересчитывать объемные параметры с геологической сетки на решетки (тоже сетки, но со значительно более крупными ячейками). В группе «Lattice» содержатся два типа объектов – решетки и кубы на этих решетках (фиг. 20). Решетки создаются через меню группы «Lattice». Как и геологические сетки, они вписываются в структурный (стратиграфический) каркас модели, но число ячеек по каждому измерению берется в несколько раз меньшее. Заполнение ячеек решетки осредненными значениями из геологического куба того или иного параметра производится при помощи диалогового окна «Ремасштабирование». Для вызова последнего необходимо выбрать в меню решетки команду «Upscale».

Фиг. 20. Объекты группы «Lattice» (гидродинамические решетки и кубы) и меню решетки

31

Объекты типа «решетки» хранятся в каталоге …\CGE\default\LATTICE. Соответствующие файлы имеют расширение «ltc». Формат файлов бинарный.

Объекты типа «куб на решетке» хранятся в каталоге …\CGE\default\ LATTICE\CUBE. Соответствующие файлы имеют расширение «lcb». Формат файлов бинарный.

1.11. Группа «Translate» (трансляторы)

Трансляторы представляют собой альтернативные прямоугольные системы координат x, y. Трансляторы ассоциируются с объектами типа «регион». Первоначально на исходном плане X,Y при помощи мыши позиционируется прямоугольный регион, причем названное позиционирование может включать поворот на заданный угол (то есть стороны региона могут не быть параллельными X и Y). Перед созданием ассоциированного транслятора регион желательно заблокировать. Для создания транслятора в меню региона выбирается команда «Задать транслятор». После этого транслятор с указанным именем («kom») появляется в группе «Translate».

Пункты меню транслятора (фиг. 21) позволяют открывать окна «План», «Пространство XYZ» и другие, отображающие объекты модели во вновь созданной системе координат.

Фиг. 21. Меню транслятора «kom»

Используя транслятор, с моделью можно работать в системе координат, ассоциированной, например, с планом 3D сейсмики. Пока трансляторы применяются редко. Транслятор «kom» создан исключительно для примера.

На этом мы завершаем знакомство с элементами модели в окне «Список объектов». Мы рассмотрели не все группы и не все типы объектов, но того, что мы рассмотрели, достаточно, чтобы двинуться дальше. Оставшиеся группы мы будем рассматривать по мере необходимости. Можно также найти описание оставшихся групп в документации.

32

2. Работа в окне «Model XYZ»

В предыдущей главе мы рассмотрели, из каких элементов состоит геологическая модель, и как эти элементы представляются в окне «Список объектов». В последующих трех главах мы будем знакомиться с тем, как элементы модели выглядят и как с ними можно работать в окнах визуализации трех основных типов - в окне 3D, в окне 2D и в окне «Профиль». Данная глава будет первой из названных трех, и в ней мы будем изучать представление элементов модели в окне 3D «Model XYZ».

После загрузки модели (фиг. 8) нажатием на кнопку «Панель инструментов» откроем кнопочную панель «Инструменты DV». При помощи команды главного меню «Пространство\Model XYZ» откроем окно с изображением контура куба. Мы получим состояние, показанное на фиг. 22.

Фиг. 22. Элементы интерфейса, позволяющие работать в окне «Model XYZ»

2.1. Действия с контуром куба

Обратим внимание на исходные границы куба по X, Y и Z. Эти границы закреплены в dvgфайле в разделе “Space”. DV-Geo позволяет работать с кубом только внутри этих границ. «Некруглые» предельные значения границ по X,Y обусловлены тем, что они обычно устанавливаются процедурой загрузки скважин по координатам последних (пользователь указывает только ширину окружающей полосы). Если исходные границы куба недостаточно широки, можно переустановить их в dvg-файле.

33

Рассмотрим, какие манипуляции возможны в окне «Model XYZ» с одним только контуром куба. Самое простое действие – используя «движки», расположенные вдоль верхней и правой границ окна, мы можем изменять размер изображения куба относительно размера окна.

Далее, можно вращать куб, установив режим того или иного вращения нажатием кнопки с соответствующей пиктограммой на панели главного окна. Вращение производится перемещением мыши с нажатой левой кнопкой по рабочей области окна «Model XYZ».

Можно сузить границы видимого куба, то есть выбрать его некоторый подкуб. Для этого необходимо нажать кнопку «Выбрать подпространство». После нажатия этой кнопки в кубе появляется селектор объема и связанное с ним диалоговое окно. Границы подкуба можно установить или перемещая границы селектора объема (удерживая их нажатой левой кнопкой мыши), или через диалоговое окно. Обратите внимание, при работе с селектором объема курсор мыши приобретает две особые формы.

Рядом с кнопкой «Выбрать подпространство» находится кнопка, позволяющая изменять пропорции куба по X,Y, Z. Желаемые изменения производятся перемещением мыши с нажатой левой кнопкой по рабочей области окна «Model XYZ». Кнопка «Отмена» приводит куб к исходным пропорциям, ракурсу и размеру в окне.

2.2. Работа мышью в окнах визуализации модели

Мы только что познакомились, как действует левая кнопка мыши в окне «Model XYZ». Результат этих действий зависит от того, какой режим интерпретации движений мыши установлен кнопками на панели главного окна и панели «Инструменты DV». Что еще можно делать мышью в окне «Model XYZ»?

Сделаем важное замечание. В системе «DV-Geo» каждое окно с трехмерным или двухмерным изображением модели (окно 3D или 2D) является двухсторонним элементом интерфейса. То есть, пользователь может не только видеть в этом окне изображения тех или иных объектов (скважин, отметок на скважинах, каротажных кривых, поверхностей, слайсов, селекторов и т.д.), но и воздействовать на видимые объекты манипулятором «мышь». Если объекты, видимые в окне 3D, проецируются друг на друга, действия мыши касаются того объекта, который к наблюдателю ближе.

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

Захват указанных мышью объектов возможен не только в окнах 3D и 2D, но и в окне «Список объектов», в окнах отображения цветов, палитр, шкал и т.д. Захватив объект, мы можем «перетащить» его в другое окно («положить» в группу, «положить» на другой объект).

34

Момент освобождения объекта от захвата называют также «бросанием» объекта (в окно, в куб, на слайс, на план, на поверхность, на линию и т.д.).

Пока в окне «Model XYZ», как нам кажется, ничего нет. Это не так. Указывая на произвольную точку внутри окна, мы указываем на объект «Пространство». Мы можем его захватить и вызвать его меню.

2.3. Отображение скважин и поверхностей

Бросим в куб поверхности «1» и «2» (то есть, захватим их пиктограммы в окне «Список объектов» и перетащим их в окно «Model XYZ»). Сделаем то же самое в отношении группы «Well». Перетаскивать поверхности и скважины из «Списка объектов» можно как группами, так и поштучно.

Результат бросания в куб поверхностей и скважин показан на фиг. 23 (слева). Скважины видны как точки на нижней грани куба. Для отображения траектории одной скважины необходимо дважды щелкнуть на соответствующей точке левой кнопкой мыши. Если при этом на панели «Инструменты DV» будет нажата кнопка «Группировка скважин», отобразятся траектории всех скважин (фиг. 23. справа). Скрываются траектории скважин повторным выполнением тех же действий.

Фиг. 23. Вид куба после помещения в него поверхностей и скважин. Показано два режима в части видимости траекторий скважин.

Реакция скважин на действия мышью в окне 3D является не совсем стандартной. Меню скважины не имеют, а на соответствующий вызов реагируют (изменением видимости траекторий) только точки пересечения скважин с нижней гранью куба. Нажатие правой кнопки мыши с указанием на скважину приводит к отображению имени скважины и координат указанной точки на траектории.

35

Поверхности в кубе отзываются на нажатия мыши стандартным образом (захватываются мышью и отображают свои меню). Для удаления из куба поверхности или скважины достаточно положить на нее букву «D» из правого нижнего угла окна «Model XYZ». Буква «D» захватывается правой кнопкой мыши. Удаление всех скважин происходит, если в момент удаления нажата кнопка «Группировка скважин».

2.4. Построение слайсов, отображение кубов параметров

Основным способом отображения рассчитанных в кубе объемных параметров (литологии, пористости и т.п.) являются слайсы. В DV-Geo можно построить слайсы (сечения) многих разных типов. Большинство из них строится при помощи кнопок на панели «Инструменты DV». Чаще всего используются две соседние кнопки, показанные стрелками на фиг. 22.

Первая из них позволяет создать так называемый «регулярный» слайс, параллельный грани куба. Нажимаем ее и указываем левой кнопкой мыши точку на одной из видимых граней куба – появляется контур регулярного слайса. Соседняя кнопка позволяет создать «произвольный» слайс, параллельный ребру куба. Нажимаем вторую кнопку и прочерчиваем мышью (при нажатой левой кнопке) отрезок на грани куба – появляется рамка произвольного слайса (фиг. 24 слева).

Фиг. 24. Три регулярных и один произвольный слайс. На правой картине на слайсы положен куб «lito».

Слайсы можно перемещать в пределах куба, захватывая их левой кнопкой мыши, но только если на панели «Инструменты DV» нажата кнопка «Перемещение» (фиг. 22). Для того чтобы слайс отображал сечение параметра, необходимо положить на него куб параметра. Раскроем в окне «Список объектов» группу «Pack». Захватим мышью куб литологии («lito») и положим его на слайс в кубе. Результат показан на фиг. 24 справа. На произвольном (вертикальном) и регулярном (горизонтальном) слайсах отображается один и тот же параметр «литология». Но на разных слайсах могут отображаться разные параметры.

Заметим следующее. Если Вы точно воспроизводили описанные действия, Вы не должны получить результат, показанный на фиг. 24. Дело в том, что область расчета литологии занимает крайне малую часть исходного объема. Для того, чтобы получить показанную

36

картину, необходимо в исходном кубе выделить фрагмент, границы которого указаны на фиг. 24 в диалоговом окне. Процедуру выделения фрагмента мы рассматривали выше.

Еще одно различие полученной Вами картины и показанной на фиг.24 справа может следовать из-за того, что у Вас с параметром «lito» может быть ассоциирована другая цветовая палитра. Палитру отображения параметра легко поменять. Для этого необходимо командой из меню главного окна «Данные\Палитры» открыть окно палитр, захватить правой кнопкой мыши любую палитру и положить ее на слайс в кубе. Палитра отображения параметра изменится. Если мы укажем мышью на окрашенную часть слайса и нажмем правую кнопку, то мы увидим, какое значение имеет параметр «lito» в той или иной точке. Параметр «lito» в нашем кубе равняется или нулю, или единице. Поэтому он при любой палитре отображается на слайсе лишь двумя цветами.

Мы видим, что реакция слайса на нажатие правой кнопки мыши является стандартной. Двойное нажатие левой кнопки мыши при указании на слайс также вызывает стандартный ответ - появление меню слайса. Состав меню зависит от типа слайса и, в некоторых случаях, от его ориентации (вертикальной или горизонтальной). Удалить слайс можно или через его меню, или положив на него букву «D» из правого нижнего угла окна. Нажатие правой кнопки мыши с указанием на слайс вне области отображения параметра ведет к захвату слайса (но только для регулярного слайса).

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

2.5. Отображение объемных тел

В окне «Model XYZ» можно показать в виде объемного тела пространственную область, в которой значения из куба параметра лежат в заданном интервале. На фиг. 25 показана пространственная область резервуара с повышенными значениями пористости. Выборку тела из исходного объема можно делать по значениям нескольких параметров. Читайте об этом «Руководство пользователя», том 5, главу 5.

37

Фиг. 25. Область резервуара с пористостью выше 22 %

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

В окне «Model XYZ» можно показать также гидродинамическую решетку и регулярные слайсы кубов на решетке (фиг. 26). Отображение решетки и слайсов производится при помощи кнопочной панели «Инструментарий гидродинамической сетки». Работа с названной панелью описывается далее в главе 6 Части II.

Фиг. 26. Гидродинамическая решетка “GRD” и регулярные слайсы куба “por” на решетке

В заключение описания окна «Model XYZ» упомянем, что в нем можно показать не только скважины, но и каротажные кривые. Однако в данном типе окна эта возможность развита слабо и ею не пользуются. Отображение данных ГИС в пространстве реализовано в новом модуле 3D визуализации (см. «Руководство пользователя», том 8). Другие ассоциированные со скважинами данные (корреляционные маркеры, разбивки РГИ и т.д.) в рассматриваемом окне показать нельзя.

38

3. Работа в окне «План горизонтального слайса»

Вернемся к состоянию модели, показанному на фиг. 24 (справа) Активируем в кубе меню горизонтального регулярного слайса, как это показано на рисунке ниже. Выберем в этом меню пункт «План». Откроется окно, ассоциированное с выбранным слайсом (фиг. 27).

Обратите внимание, во вновь открытом окне отображается то же самое сечение куба «lito», что и на слайсе в кубе. При перемещении слайса внутри куба будет меняться и картина в ассоциированном окне.

Фиг. 27. Окно регулярного горизонтального слайса

Окно «План горизонтального слайса» относится ко второму важному типу окон отображения модели – к типу окон 2D. Вообще в DV-Geo есть много разновидностей окон 2D. Но при этом именно окно горизонтального слайса является по функциональности одним из главных.

Обратите внимание ниже на размер его меню (фиг. 29). Почти таким же (но с несколько меньшей функциональностью) является окно «План XY», которое открывается командой главного меню «Пространство\План XY». Рассмотрим, какое представление элементы модели имеют на плане горизонтального слайса.

3.1. Настройка размера шрифта на плане XY

В окне «План горизонтального слайса» отображаются координатные оси X,Y. В начальный момент на Вашем экране шрифт разметки координатных осей может быть слишком крупным или, наоборот, слишком мелким (как фиг. 27). Координатные оси отзываются на воздействие мышью стандартным образом, то есть могут показать свое меню. Через это меню можно изменить размер шрифта. Но не стоит торопиться и действовать именно таким образом, причина может быть в другом. А именно, возможно, что размер Вашего экрана существенно отличается от размера экрана, на котором модель визуализировалась ранее. Вследствие этого обстоятельства неподходящими могут оказаться очень многие установки. Для проверки соответствия текущих установок размеру экрана необходимо выбрать пункт главного меню «Файл\Твердая копия\Новая». Откроется окно формы, показанное на фиг. 28. В этом окне необходимо нажать кнопку «Размер экрана». Откроется одноименное диалоговое окно. Необходимо сравнить числа в полях этого окна с фактическими размерами рабочей области

39

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

Фиг. 28. Установка размеров окна

3.2. Работа мышью на плане XY

Захватим мышью в «Списке объектов» группу скважин Well и бросим ее на план (горизонтального) слайса. Получим картину, показанную на фиг. 29. Таким же образом на план слайса можно бросить и многие другие объекты – кубы параметров, карты, изолинии, объекты типа Area, гидродинамические кубы и т.д.

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

40

Разобраться с тем, какие объекты присутствуют на плане, помогает правая кнопка мыши. При ее нажатии с указанием на тот или иной объект или высвечивается краткое информационное сообщение об объекте (имя для скважины, координата для оси, значение параметра для сечения), или производится захват объекта (для плана). Это позволяет понять, меню какого объекта всплывает в момент двойного щелчка левой кнопкой мыши. На фиг. 29 представлены меню плана, сечения параметра и скважины, и указаны точки, в которых производились вызовы соответствующих меню. Ранее было показано меню координатных осей (фиг. 27).

На плане XY реализована стандартная функция зуммирования. Условием ее применения является отжатое положение кнопки «Перемещение» на кнопочной панели «Инструменты DV». Вычетрим мышью с нажатой левой кнопкой прямоугольник на плане XY (движеним мыши вправо-вниз). При отпускании левой кнопки мыши произойдет расширение выделенного прямоугольника до размеров рабочей области окна. То же действие, при котором прямоугольник рисуется обратным перемещением влево-вверх, производит противоположный эффект - текущая рабочая область окна умещается в контур выделенного прямоугольника (пока не достигаются максимальные значения координат по внешним границам). В режиме «Перемещение» движение мыши с нажатой левой кнопкой производит сдвиг области зуммирования на плане X,Y (если зуммирование имеет место).

3.3. Отображение скважин

Мы уже визуализировали на плане расположение скважин (фиг. 29). Возможно, у Вас скважины отобразятся без номеров. В этом случае надо вызвать меню плана (показано на фиг. 29 слева) и выбрать в нем пункт «Имена скважин\Имена всегда». Изменить размер шрифта номеров скважин позволяет соседний пункт того же меню. Величину отметки (кружка) скважины можно изменить через меню скважины. Пункт «Свойства» меню скважины открывает диалоговое окно «Свойства скважин». Размер кружка определяется значением в первом из двух полей против имени «Диаметр скважины».

Какие координаты скважин мы видим на плане горизонтального слайса, устьев или забоев? Если мы хотим видеть на плане точки пересечений траекторий скважин с плоскостью слайса, мы должны в меню плана активировать пункт «План по глубинам». Названный пункт меню отмечается галочкой, после чего к отметке каждой невертикальной скважины добавляется некоторый отрезок (фиг. 30). Отметка (кружок) имеет координаты пересечения скважины с плоскостью слайса, а противоположный конец отрезка имеет координаты забоя скважины. При перемещении слайса в кубе по вертикали, точки его пересечения со скважинами на плане тоже перемещаются, но только пока глубина слайса меньше глубины забоя. Ниже забоя каждая скважина условно продолжается вниз вертикально. В случае отмены опции «План по глубинам» отметки скважин будут иметь координаты забоев.

41

Фиг. 30. Отображение на плане невертикальных скважин

Часто в пределах целевого горизонта используется вертикальная аппроксимация траекторий скважин. При этом в качестве координат X,Y скважин на плане берутся координаты пластопересечений. Оговорка «в пределах целевого горизонта» позволяет использовать вертикальную аппроксимацию даже для кустов скважин. В проекте «cge» используется именно такая аппроксимация. В этом случае опция «План по глубинам» ничего не меняет. Иллюстрация, показанная на фиг. 30, касается другого проекта.

Выше мы упоминали, что скважина на плане (в отличие от скважины в окне 3D) имеет свое меню. Обратите внимание, последнее отличается от меню скважины в окне «Список объектов». Но ряд возможностей дублируется. Так, через пункт меню «Список объектов» скважины на плане можно получить точно такой же доступ ко всем данным соответствующей скважины, как и через полный «Список объектов». Пункт «Открыть» в меню скважины на плане позволяет визуализировать скважинные данные в окне типа «Профиль» (окно этого типа будет описано дальше). Пункт «Блокировать» позволяет исключить данные скважины из всех последующих расчетов (интерполяции поверхностей, интерполяции кубов и т.п.).

3.4. Использование фильтра скважин

Важной функций окна «План (горизонтального) слайса» является фильтрация скважин по указанному признаку. Вызов диалогового окна «Фильтр» производится из меню плана XY (фиг. 31). На иллюстрации показана наиболее простая фильтрация по признаку наличия на скважинах отметок корреляционного репера. Для примера выбран репер «bot73c». Названный репер переносится из списка объектов в поле фильтра «Первый» (имеется в виду репер, можно указать и второй), после чего на панели фильтра нажимается кнопка «Использовать фильтр». В результате на плане XY остаются только скважины, на которых имеются отметки указанного репера (25 скважин из 29). После нажатия кнопки «Инвертировать» на плане отобразятся 4 скважины, на которых отметки данного репера отсутствуют.

42

Фиг. 31. Фильтрация скважин на плане XY

Далее мы познакомимся с так называемыми селекторами скважин на плане XY и увидим, что селекторы отбирают только скважины, прошедшие заданный фильтр. Фильтрация скважин по наличию отметок корреляционного репера помогает при выполнении ручной корреляции скважин по данным ГИС.

Фильтр позволяет не только выделять скважины по тому или иному признаку, но и блокировать выделенные скважины. Для этого необходимо нажать кнопку «Блокировать». Напомним, что блокированные скважины не принимаются во внимание при выполнении расчетных процедур.

3.5. Отображение регионов, поверхностей, карт

Важной функций окна «План (горизонтального) слайса» является редактирование регионов для последующего расчета поверхностей и кубов. Мы уже видели, что в модели «CGE» используется только один регион «1». Если захватить его мышью в окне «Список объектов» и бросить на план XY, то отобразится его контур (фиг. 32). Контур региона на имеет меню, которое также показано на фиг. 32. Через это меню можно вызвать окно с параметрами региона и посмотреть размер ячеек его сетки. Но, поскольку на регионе уже созданы поверхности, изменить размер ячеек его сетки уже нельзя.

Новый регион можно создать через меню группы «Region» окна «Список объектов». При создании он получает размер, близкий к максимально возможному. Изменить размеры региона (пока на нем не созданы поверхности) можно, захватывая его границы на плане XY левой кнопкой мыши. Захват границ (как и вызов меню региона) производится с их внутренней стороны. Регион можно вращать на плане, захватывая его за углы.

43

Фиг. 32. Вид региона «1», меню региона и диалоговое окно региона «Параметры»

Что важно учитывать при выборе размера ячеек сетки нового региона? С одной стороны, необходимо помнить, что сеть на регионе будет определять расчет не только поверхностей, но и кубов. Куб 200х200х200 (8 миллионов ячеек) занимает более 32 Мб. Но, с другой стороны, сеть не должна допускать попадания двух скважин в одну ячейку, поскольку в таком случае одну из двух скважин приходится исключать (блокировать). В диалоговом окне «Регион» можно задать режим отображения сетки, а также режимы «Одна ячейка» и «Соседние ячейки». В режиме «Одна ячейка» выделяются ячейки, в которые попадают две скважины. В режиме «Соседние ячейки» выделяются ячейки со скважинами, соседствующие друг с другом.

Фиг. 33. Регион «1» и поверхность «1» модели «CGE»

На фиг. 33 (слева) регион «1» показан в режиме отображения сетки. Можно видеть, что контур региона определен таким образом, чтобы с небольшим запасом охватить все скважины. Размер ячейки по X и Y выбран равным 50 м. Размер сетки равен 101х113 ячеек.

44

О том, как на регионе создается поверхность, мы говорили при рассмотрении региона и поверхности в окне «Список объектов». После расчета поверхности параметры региона блокируются, то есть изменять их уже нельзя. Бросим в окно «План слайса» поверхность «1». Появится меню, в котором необходимо выбрать или опцию «Карта», или опцию «Пересечение». Выбор опции «Карта» даст вид поверхности, показанный на фиг. 33 в центре. Альтернативная опция «Пересечение» дает отображение линии пересечения поверхности и слайса. На той же иллюстрации справа показано окно с характеристиками сетки поверхности «1», отображаемое через меню поверхности в «Списке объектов» (пункт «Информация»). Число значений поверхности 102х114, поскольку, как мы говорили, они ассоциированы не с центрами ячеек, а с узлами сетки.

Как, имея карту поверхности, рассчитать изолинии? В DV-Geo есть две возможности построения изолиний. Первая возможность ручная, интерактивная. Командой главного меню «Данные\Цвета» открываем окно с палитрой цветов. Правой кнопкой мыши захватываем любой цвет (это будет цвет изолинии) и переносим его на карту на плане слайса, в точку с некоторым значением (это будет уровень изолинии). В результате этих действий на карте появляется одна изолиния (фиг. 34). Далее выбираем в меню поверхности на плане пункт «Шкала» и открываем шкалу с отметкой рассчитанной изолинии (фиг. 34 справа). Эту отметку можно двигать мышью, изолиния на карте будет изменяться синхронно. Двойной щелчок левой кнопки мыши на той же отметке открывает окно «Установить уровень», в котором можно задать начальное значение и шаг для расчета набора изолиний. Особенность этих изолиний в том, что их нельзя сохранить, и они пропадают при закрытии карты.

Фиг. 34. Построение изолиний посредством выделения цветом уровней на карте

Вторая возможность расчета изолиний – использовать пункт «Создать карту в изолиниях» в меню поверхности в «Списке объектов». Изолинии, рассчитанные таким образом, сохраняются как все другие объекты модели. Смотрите также на данную тему «Руководство пользователя», том 2, раздел «Построение карт в изолиниях».

3.6. Отображение кубов параметров. Седиментационный слайс

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

45

положить и другой куб параметра, то есть кубы на слайсе в пространстве 3D и в ассоциированном окне могут быть разными. Проблема в том, что сечение куба параметра горизонтальной плоскостью не является информативным, поскольку горизонтальная плоскость сечет стратиграфические горизонты. Сделать горизонтальный слайс информативным можно двумя путями. Первый путь – использовать не строго горизонтальный, а так называемый «седиментационный» слайс, идущий параллельно слоям стратиграфической сетки. Второй путь – выполнить палеореконструкцию модели, то есть превратить структурные поверхности модели в горизонтальные плоскости.

Построение седиментационного слайса показано на фиг. 35. После нажатия кнопки «Седиментационный анализ» на кнопочной панели главного окна DV-Geo открывается окно «Шкала». На противоположные края этой шкалы необходимо положить поверхности «1» и «2», перенеся их из «Списка объектов». Названия поверхностей отобразятся на концах шкалы. После этого красная отметка на шкале будет определять поверхность, которая делит интервал между «1» и «2» в заданной пропорции. Отметку на шкале можно перемещать мышью. Если захватить эту отметку правой кнопкой мыши и бросить в окно 3D, в кубе отобразится соответствующая поверхность. Двойной щелчок на отметке левой кнопкой мыши приведет к открытию окна «План для седиментационного анализа». Если положить в это окно куб «lito» (или любой другой), отобразится соответствующее сечение. Можно положить куб «lito» и на поверхность седиментационного слайса в пространстве 3D. Сечение куба седиментационным слайсом отображает черты обстановки осадконакопления, то есть является информативным.

Фиг. 35. Седиментационный слайс куба литологии

46

По своим характеристикам окно седиментационного слайса аналогично окну горизонтального слайса. Для лучшего сходства картины в окне 2D с картиной на поверхности слайса в кубе мы поменяли местами оси на плане XY (через меню этих осей) и инвертировали направление оси X.

Покажем теперь, как ту же картину можно получить с помощью палеореконструкции модели. Идея палеореконструкции исходит из того, что стратиграфические горизонты в момент своего формирования являлись горизонтальными. Следовательно, если мы выровняем два соседних горизонта, мы восстановим черты древней обстановки осадконакопления. Горизонтальный слайс между выровненными поверхностями станет при этом полным аналогом седиментационного. На фиг. 36 показан горизонтальный слайс параметра «lito» в условиях, когда обе помещенные в куб поверхности «1» и «2» превращены в горизонтальные плоскости и сделаны «прозрачными». Как мы и ожидали, картины на горизонтальном (в условиях палеореконструкции) и седиментационном слайсе совпадают. Режим палеореконструкции модели по той или иной поверхности задается через меню этой поверхности в окне 3D. Степень прозрачности поверхности в кубе устанавливается через то же меню.

Фиг. 36. Горизонтальный слайс куба литологии в условиях палеореконструкции модели

3.7. Отображение объектов типа «Area»

Для отображения на плане XY (например, в окне горизонтального слайса) объекта типа «Area», необходимо захватить его мышью в списке объектов и перенести на названный план. Объект типа «Area» изображается на плане как линия, замкнутая или не замкнутая. На фиг. 37 показано изображение на плане объекта «raz5_2_1» (красная линия). Это контур разлома. Вместо меню линия на плане открывает диалоговое окно, с помощью которого можно изменить ее толщину, цвет, стиль, отобразить точки и т.д. (фиг. 37 справа). Если отметить в этом окне на вкладке «Общие» пункт «Редактировать», появится дополнительная кнопочная панель, которая позволит изменить саму линию.

47

Фиг. 37. Отображение на плане XY объекта типа «Area»

Аналогичным образом можно отобразить и другие объекты типа «Area». Об их назначении мы говорили при рассмотрении «Списка объектов».

Объекты «Area» разделяются на ряд подтипов – замкнутые и незамкнутые произвольные линии, линии через скважины, овальные и прямоугольные селекторы скважин. Для создания на плане нового объекта типа «Area» необходимо нажать одну из шести кнопок на панели «Инструменты DV» (фиг. 38). Различия между кнопками понятны из нанесенных на них пиктограмм и поясняющих надписей на рисунке. В установленном режиме на плане XY мышью с нажатой левой кнопкой вычерчивается или соответствующая линия, или селектор. После того, как построен первый отрезок линии, можно захватить мышью одну из его внутренних точек и превратит отрезок в замкнутый или незамкнутый контур. Построенный селектор также можно захватывать левой кнопкой мыши и, в зависимости от точки захвата, перемещать, вращать, сжимать и растягивать.

Для сохранения созданного объекта (линии, селектора) необходимо захватить его правой кнопкой мыши и положить на название группы «Area» (или на название одной из вложенных подгрупп) в «Списке объектов». В результате этого действия появляется диалоговое окно «Новая область», в которое необходимо ввести имя области. После введения имени пиктограмма объекта появляется в группе «Area» (или в соответствующей подгруппе).

Прямоугольные селекторы и линии, проходящие через скважины (в отличие от просто линий), имеют свои меню. В этих меню есть пункт «Профиль». Данный пункт позволяет открыть окно типа «Профиль», в котором отображаются вертикальные разрезы по включенным скважинам. Для овального селектора окно «Профиль» открывается пунктом меню «Открыть». Окно типа «Профиль», является в DV-Geo одним из главных рабочих окон. Его рассмотрению будет отведена целая глава.

48

Фиг. 38. Кнопки режимов для создания новых объектов типа «Area»

49

4. Работа в окне «План вертикального слайса»

Вернемся к состоянию модели, показанному на фиг. 24 (слева) Активируем в кубе меню вертикального регулярного слайса, как это показано на рисунке ниже. Выберем в этом меню пункт «План» (но не «План слайса»). Откроется окно, ассоциированное с выбранным слайсом (фиг. 39). Во вновь открытом окне будет отображено то же самое сечение куба «lito», что и на слайсе в кубе. При перемещении слайса внутри куба будет меняться и картина в ассоциированном окне.

Фиг. 39. Окно регулярного вертикального слайса

4.1.Использование окна вертикального слайса для переустановки границ куба в окне

«Model XYZ»

Окно регулярного вертикального слайса удобно использовать для выделения отображаемого фрагмента куба не только в нем самом, но и в окне «Model XYZ». Действительно, стандартная функция зуммирования в окне «План слайса» (см. раздел 3.2) позволяет быстро и точно увеличить область окна, представляющую интерес. После этого границы по осям Y и Z на плане можно захватить мышью и перенести в окно «Model XYZ». Захват осей производится указанием на одну из них левой кнопкой мыши (так, чтобы появилось значение координаты) и последующим дополнительным нажатием правой кнопки. При этом рядом с курсором мыши появляется рамка с текстом «Оси плана», которую необходимо перетащить в окно «Model XYZ». Результат показан на фиг. 40. Границы куба по оси X при этом не изменяются.

50

Фиг. 40. Изменение границ куба при помощи плана вертикального слайса

51

5. Работа в окне «Профиль»

Создадим в окне «План XY» прямоугольный селектор скважин. Напомним, как это делается.

Сначала нажатием кнопки на панели «Инструменты DV» (см. также фиг. 38), установим соответствующий режим интерпретации действий мыши. Затем вычертим селектор мышью (с нажатой левой кнопкой) на плане XY. Потом вызовем меню созданного селектора и в этом меню выберем пункт «Профиль». Откроется окно «Профиль», показанное на фиг. 41.

Фиг. 41. Окно «Профиль», ассоциированное с прямоугольным селектором скважин

Точно такое же окно «Профиль» можно открыть из меню линии, проходящей через скважины, и из меню овального селектора. Окно «Профиль» можно открыть и непосредственно через меню главного окна DV-Geo, выбрав пункт «Пространство/Профиль». Рассмотрим, что собой представляет окно «Профиль» и какие элементы модели можно в нем показать.

5.1. Общее описание окна «Профиль»

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

Какие именно скважины и в каком порядке представлены в окне «Профиль»? Для окна, ассоциированного с линией профиля через скважины, ответ очевиден. Если это окно, ассоциированное с прямоугольным селектором, то в нем представлены скважины, попадающие в границы селектора, и в том порядке, в каком они проецируются на более

52

длинную сторону селектора. При перемещении селектора состав и порядок отображения скважин синхронно меняется. Если окно «Профиль» открыто через меню главного окна, скважины переносятся в него мышью из списка объектов.

Окно «Профиль» имеет четыре варианта отображения фреймов скважин – «Полный», «Простой», «Проекция», «Профиль». На фиг. 41 показан вариант «Профиль». Вариант можно переустановить через меню шкалы глубин или через меню любой скважины (названные меню и примерные точки их вызова показаны на фиг. 41). Опишем коротко эти варианты. «Полный» - фреймы скважин имеют одинаковую ширину, каждая скважина имеет свою шкалу глубин и может масштабироваться по этой шкале независимо от остальных. «Простой» - шкала глубин у скважин общая, но фреймы всех скважин также имеют одинаковую ширину. «Профиль» - скважины внутри селектора условно соединяются ломаной линией по нарастанию координаты X, и ширина фрейма каждой скважины (кроме последней) пропорциональна длине исходящего отрезка ломаной. «Проекция» - ширина фреймов определяется проекциями скважин на сторону селектора X. В двух последних случаях отображение данных скважины может выходить на соседние фреймы.

Поскольку скважинные данные имеют два вида координатной привязки – первичную кабельную длину и расчетные пространственные координаты X,Y, Z, их отображение в окне «Профиль» возможно и в координатах кабельной длины, и в координатах абсолютной глубины. В верхней части шкалы глубин стоит буквенное обозначение текущего режима - D (абсолютная глубина) или L (кабельная длина). Есть также два дополнительных режима - DL (абсолютная глубина плюс отображение графика удлинения кабеля) и T. В последнем случае вместо глубины отображается время (даты). Буквенное обозначение имеет меню (фиг. 41), через которое можно выбрать нужный режим.

Масштабирование по шкале глубин можно производить разными действиями. Можно захватывать осевую линию шкалы левой кнопкой мыши и смещать ее вверх-вниз. Если при этом нажата кнопка «Перемещение» на панели «Инструменты DV» (фиг. 22), шкала будет прокручиваться вверх-вниз. Если же названная кнопка отжата, шкала будет растягиваться или сжиматься при фиксированной верхней отметке. Другая возможность масштабирования по глубине – вычертить мышью прямоугольник, охватывающий фрагмент шкалы (фиг. 41). Если прямоугольник вычерчивается движением «слева направо и вниз», выделенный фрагмент будет растянут на всю шкалу. Если же прямоугольник вычерчивается движением «слева направо и вверх», текущий интервал шкалы будет вмещен в выделенный фрагмент.

Прямоугольный селектор скважин можно перемещать не только мышью на плане XY, но также используя кнопки перемещения селектора на кнопочной панели окна «Профиль» (фиг. 41). Перемещение селектора производится в направлении, перпендикулярном его длинной стороне (вперед, назад и автоматически в сторону последнего единичного шага). Один шаг перемещения селектора составляет 1/3 его ширины.

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

5.2. Отображение и редактирование данных ГИС

Окно «Профиль» насыщено сервисами, касающимися визуализации и редактирования каротажных кривых. Для того, чтобы на всех скважинах в окне «Профиль» визуализировать

53

каротажные кривые одного метода ГИС, необходимо захватить каротажную кривую этого метода в списке объектов (в группе любой скважины) и перенести ее на шкалу глубин окна «Профиль». После этого каротажные кривые данного метода отобразятся на всех скважинах, на которых они имеются (фиг. 42). Если название метода ГИС бросить не на шкалу, а на фреймы скважин, соответствующая кривая отобразится только на «своей» скважине (на которой она была захвачена).

Фиг. 42. Отображение и редактирование каротажных кривых в окне «Профиль»

Список объектов, имеющихся на той или иной скважине, можно открыть через меню скважины в окне «Профиль». При помещении каротажных кривых в окно «Профиль» не имеет значения, где они были захвачены – в общем окне «Список объектов», или в списке объектов одной скважины.

Каждая каротажная кривая в окне «Профиль» имеет свое меню (также показано на фиг. 42). Назовем лишь некоторые его функции. Пункт «Удалить» касается только изображения кривых данного типа в окне, но не удаления самих кривых из модели. Через пункт «Свойства» можно изменить толщину кривых. Через пункт «Средства» пользователь может вставить в это меню собственную расчетную процедуру, и она будет получать имя той самой кривой, которой принадлежит меню.

Для изменения цвета каротажной кривой необходимо через меню главного окна DV-Geo «Данные/Цвета» открыть окно «Цвета», захватить левой кнопкой мыши желаемый цвет и положить его на выбранную каротажную кривую. Если при этом на панели «Инструменты DV» будет нажата кнопка «Группирование однотипных ГИС», заданный цвет приобретет не только выбранная кривая, но и все каротажные кривые того же типа. Забегая вперед, скажем,

54

что таким же образом в окне «Профиль» можно изменять цвет разбивок РГИ, отметок на скважинах и линий пересечения поверхностей.

Для изменения амплитуды отображения каротажной кривой необходимо установить режим «Изменение амплитуды» на кнопочной панели «Инструменты DV» (фиг. 42). После этого необходимо захватить левой кнопкой мыши каротажную кривую и изменить ее амплитуду. Если в этот момент нажата кнопка «Группирование однотипных ГИС» (или если скважины отображаются в режимах «Профиль», «Проекция»), изменение амплитуды будет касаться всех однотипных каротажных кривых. При нажатой кнопке «Группирование скважин» изменение будет касаться каротажных кривых всех типов на всех скважинах. Кнопка «Группирование ГИС одной скважины» распространяет изменение на все типы кривых одной скважины. Названные действия не меняют значения каротажных кривых.

Для смещения каротажной кривой влево или вправо необходимо установить на той же кнопочной панели режим «Изменение выноса». Все остальные действия производятся аналогично изменению амплитуд.

На кнопочной панели окна «Профиль» имеется кнопка, устанавливающая режим ручного редактирования кривых каротажа (фиг. 42). В этом режиме левой кнопкой мыши можно вручную перерисовать (исправить) любой фрагмент. Новая линия должна начинаться и заканчиваться пересечением с редактируемой кривой. Если пересечение только одно, оно интерпретируется как команда обрезать исходную линию, выше или ниже точки пересечения. Появляется окно с предложением сделать указанный выбор, а затем, при получении того или иного утвердительного ответа, происходит удаление соответствующей части кривой.

Важной редакторской функцией является увязка каротажных кривых. Данные разных методов ГИС могут быть получены как в ходе одного «прохода» (подъема на кабеле комплексного прибора), так и в ходе разных «проходов». В последнем случае может иметь место рассогласование кривых, то есть ситуация, когда одинаковые зарегистрированные глубины двух измерений на самом деле соответствуют разным фактическим глубинам.

Факт рассогласования кривых устанавливает геолог. Он же может устранить рассогласование кривых на скважине, используя режим «Увязка кривых» (соответствующая кнопка показана на фиг. 42).

55

Фиг. 43. Увязка каротажных кривых

Пример увязки кривых показан на фиг. 43. Обратите внимание, фреймы скважин показаны в режиме «Простой». Допустим, на скважине 644 мы желаем согласованно сместить кривые IK и GK (синюю и черную) относительно кривой PS (красной). Сначала скрепим кривые IK и GK между собой. Это делается при помощи горизонтальной шкалы. Нажатием правой кнопки мыши захватим в верхней части фрейма имя скважины Well 644 и бросим его на тот же фрейм. В нижней части фрейма появится горизонтальная шкала. Захватим правой кнопкой мыши горизонтальную шкалу и бросим ее на кривую GK. На шкале появится значение GK. Бросим шкалу на кривую IK, на шкале появится значение IK. После этого кривые IK и GK будут увязываться согласованно. Открепить кривые друг от друга можно повторением тех же действий. Далее необходимо указать на кривых IK и PS точки, имеющие равную фактическую глубину. Захватим левой кнопкой мыши точку на кривой IK и протянем «резиновую нить» до соответствующей ей точки на кривой PS. Теперь необходимо выровнять кривые. Вызовем меню той части «резиновой нити», которая ближе к IK, и выберем команду «Выравнивание IK». Кривые IK и GK сместятся вниз, так, что «резиновая нить» станет горизонтальной. Перед выравниванием можно протянуть и вторую резиновую нить, поставив в соответствие друг другу вторую пару точек на кривых IK и PS. В этом случае по команде «Выравнивание IK» произойдет выравнивание обеих «резиновых нитей». Команда «Сохранить проект» зафиксирует результат увязки кривых.

В полном объеме возможности редактирования кривых ГИС в окне «Профиль» описаны в «Руководстве пользователя», том 6, раздел 13.

5.3. Отображение корреляционных реперов

Мы уже знаем, что результаты корреляции скважин по данным ГИС находятся в группе «Correlation» окна «Список объектов» (см. раздел 1.5). Это четыре репера, кровли и подошвы

56

горизонтов 71с и 72с. Бросим их поочередно в окно «Профиль». Полученная картина показана на фиг. 44. Для определенности показано также положение селектора на плане XY.

Фиг. 44. Отображение корреляционных реперов в модели «CGE»

Каждый репер, показанный в окне «Профиль», имеет свое меню (также показано на фиг. 44). Действие большинства его пунктов очевидно. Отметим важное. Пункт «Удалить» касается удаления только отображений отметок репера на профиле, но не самих отметок из модели. Восстановить отображения отметок можно повторным бросанием репера в окно. Напротив, пункты «Очистить репер», «Очистить репера профиля», «Очистить все репера» действительно удаляют отметку репера, соответственно, или на одной скважине, или на всех скважинах данного профиля, или на всех скважинах модели. Если после команды «Очистить…» выполнить сохранение модели, удаление отметок репера станет необратимым.

Пункт «Выравнивание» выстраивает скважины вровень по отметкам соответствующего репера. При этом на шкале у линии репера появляется треугольная отметка уровня. Одновременно можно выровнять скважины на профиле по двум реперам. Однако, при выравнивании еще и по второму реперу, интервалы скважин между реперами будут деформированы (сжаты или растянуты). На фиг. 44 показано выравнивание скважин именно по двум реперам (обратите внимание на вид отметок выравнивания на шкале). Опция «Выравнивание» значительно облегчает корреляцию скважин при наличии опорных реперов вблизи коррелируемых границ.

Перед установкой на скважинах нового репера необходимо создать его в «Списке объектов». Это делается через меню группы «Correlation». Затем вновь созданный (но еще не проставленный) репер необходимо бросить в окно «Профиль». При этом имя репера появится в выпадающем списке на кнопочной панели названного окна. Обратите внимание, в списке присутствуют имена всех отображаемых в окне реперов (поскольку их тоже бросали в окно «Профиль»), но видимым именем является только одно. Видимое имя определяет, какой именно репер устанавливается (или корректируется) пользователем в данный момент.

57

Режим установки (корректировки) репера задается нажатием кнопки «Двигать метку» (фиг. 44). После этого для установки (корректировки) репера на скважине необходимо указать мышью на нужный уровень в середине фрейма соответствующей скважины, нажать (не отпуская) левую кнопку мыши и чуть-чуть сдвинуть мышь. Поперек фрейма появится соответствующая отметка.

Режим установки отметок имеет следующую особенность. Если в момент нажатия левой кнопки мыши мы (случайно или намеренно) попадем на отметку другого репера, произойдет автоматическое переключение имени устанавливаемого репера. Эта особенность позволяет быстро переключаться с исправления одного репера на исправление другого. Однако начинающий пользователь должен внимательно следить за этими переключениями, то есть за тем, какой именно репер он устанавливает в данный момент. Часто лучше сначала установить ту или иную отметку на свободном месте, а затем, не отпуская левую кнопку мыши, переместить ее на нужную глубину. Необходимо также отключать режим «Двигать метку» сразу после того, как в нем отпадает необходимость, иначе можно нечаянно сбить точное положение той или иной отметки.

Отметки репера и даже глубину этих отметок можно отобразить на плане XY. Для этого необходимо:

1.Отметить в меню репера в окне «Профиль» пункт «Предварительный просмотр»;

2.Настроить должным образом ассоциированную с репером цветную палитру глубин;

3.Захватить репер в «Списке объектов» и бросить его на план XY;

4.Через меню скважины на плане XY и диалоговое окно «Свойства скважин» установить оптимальный размер схемы отображения репера.

После этого отметки репера будут отображаться на плане XY посредством схем (колец), раскрашенных в соответствии с их глубиной (фиг. 45). Такое отображение позволяет видеть, на каких скважинах есть отметка репера, а на каких ее нет, и представлять положение отметок по глубине.

Фиг. 45. Отображение глубины репера на плане горизонтального слайса

58

Что значит «настроить палитру глубин»? Палитра репера (имя которого стоит в выпадающем списке) в обязательном порядке присутствует на шкале глубин окна «Профиль» (фиг. 45). На ее место можно положить любую другую палитру из списка палитр. Изначально палитра может не соответствовать ни средней глубине репера, ни диапазону его глубин . Это несоответствие необходимо устранить. Палитру можно перемещать по шкале глубин, захватывая ее левой кнопкой мыши, а также растягивать за концы или, наоборот, сжимать.

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

О других возможностях, предоставляемых пользователю в окне «Профиль», смотрите «Руководство пользователя», том 3-1 «Корреляция разрезов скважин».

5.4. Отображение данных РГИ, колонок литологии и насыщения, поверхностей и кубов параметров

Отображение и работа в окне «Профиль» с данными РГИ практически полностью аналогична работе с данными ГИС. Общие принципы сохраняются и при отображении колонок литологии и насыщения.

В окне «Профиль» можно отобразить поверхность. Для этого необходимо захватить поверхность в группе Region окна «Список объектов» и бросить ее в окно «Профиль». Результат отображения поверхности «1» на фиг. 46. Мы уже знаем, что в режиме «Профиль» мы видим в окне развертку вертикального сечения вдоль ломаной линии, проходящей через отметки скважин на плане XY внутри селектора (или вдоль линии, заданной пользователем). Соответственно, поверхность в окне «Профиль» мы видим как линию пересечения на этой развертке.

Фиг. 46. Отображение поверхности «1» в окне «Профиль»

59

По линии пересечения в окне «Профиль» очень трудно судить о рельефе поверхности. Но по ней очень легко судить о том, насколько точно поверхность проходит через исходные корреляционные отметки на скважинах. Обратите внимание – поверхность «1» точно проходит через отметки репера «top71c». Линия поверхности имеет свое меню. В диалоге «Свойства» стиль линии «L20S20» означает 2 миллиметра линия (Line), 2 миллиметра промежуток (Space).

В окне «Профиль» можно отобразить и куб параметра. Для этого необходимо захватить куб в группе Pack окна «Список объектов» и бросить его в окно «Профиль». Результат отображения куба «lito» показан на фиг. 47. Здесь, повторим, мы видим развертку вертикального сечения куба «lito» вдоль ломаной линии, проходящей через скважины. Обратите внимание - значения в кубе «lito» очень хорошо воспроизводят исходные значения колонки «Литология» в скважинах. Второе, на что необходимо обратить внимание – не имеет ли параметр «lito» каких-нибудь необъяснимых аномалий в пространстве. Никаких особенных аномалий в кубе параметра «lito» мы не видим.

Фиг. 47. Отображение куба литологии в окне «Профиль»

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

60

Фиг. 48. Отображение куба пористости в окне «Профиль»

Фиг. 49. Отображение куба насыщения в окне «Профиль»

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

61

5.5. Шаблон визуализации данных в окне «Профиль»

На фиг. 50 показана дополнительная возможность оформления окна «Профиль». В меню корреляционной отметки был выбран пункт «Установить цвет интервала», который позволил раскрасить интервалы между реперами. Ту же самую раскраску можно выполнить, захватывая тот или иной цвет в окне «Цвета» и бросая его на пространство между реперами.

Фиг. 50. Сохранение шаблона визуализации данных в группе «Patterns/Profiles»

После того, как пользователь перенес в окно «Профиль» все необходимые ему объекты и настроил параметры отображения этих объектов (интервал глубин, набор кривых, амплитуды и смещения кривых, репера, поверхности, их цвета, палитры, толщины, стили и т.д.), возникает проблема быстрого воспроизведения того же вида окна во время следующего сеанса работы с DV-Geo. Или, возможно, мы захотим открыть в том же виде окно другого профиля. Эта проблема решается посредством сохранения шаблона визуализации окна «Профиль» в списке объектов. Для хранения названных шаблонов служит группа

«Patterns/Profiles» (фиг. 50).

Первое, что необходимо сделать для сохранения шаблона окна «Профиль» – добавить шаблон в группу «Profiles». Это делается через меню названной группы. При добавлении шаблона необходимо указать его имя, в примере было указано имя «схема 1». Затем необходимо вызвать меню элемента «схема 1» и выбрать в нем пункт «Редактировать». Откроется диалоговое окно для описания шаблона, показанное на фиг. 51. Однако заполнять его поля вручную нет никакой необходимости. В нижней части этого окна есть клавиша «Установить». Достаточно захватить ее правой кнопкой мыши и перенести в то окно «Профиль», которое подлежит описанию. После этого все поля шаблона заполнятся автоматически. Последнее, что необходимо сделать для воспроизведения окна «Профиль» – определить шаблон «схема 1» (через его меню) в качестве текущего.

62

Фиг. 51. Диалоговое окно для определения шаблона окна «Профиль»

На этом мы заканчиваем рассмотрение окна «Профиль». Разумеется, мы рассмотрели только основные возможности работы в этом окне. О других возможностях смотрите «Руководство пользователя», в том 2, раздел 4-3.

Подведем итог всей Части I. Итак, что новое мы узнали?

Первое - мы узнали, из каких основных объектов состоит геологическая модель резервуара нефти и газа. Это скважины и скважинные данные (инклинометрии, ГИС, РГИ, корреляционные отметки и т.д.), линии и контуры «Area», поверхности (структурные и флюидных контактов), геологические и гидродинамические кубы. Мы также узнали, какие вспомогательные объекты используются при построении и отображении модели. Это линии профилей, слайсы, селекторы, разного рода таблицы и шаблоны.

Второе - мы узнали, как элементы модели и вспомогательные объекты отображаются в окнах четырех основных типов - «Список объектов», «Model XYZ», «План XY», «Профиль». То есть, мы достаточно подробно познакомились с интерфейсом системы DV-Geo. Мы также получили представление о взаимосвязи окон отображения разных типов и о методике их совместного использования.

Третье - мы познакомились с процедурной частью системы и узнали, как процедуры встраиваются в группу «Tools» окна «Список объектов» и в меню объектов модели.

Чего мы пока не знаем? Конечно, мы пока не знаем многих частных деталей, относящихся к первому, второму и третьему. Но это не важно. Даже специалисты, долгое время работающие с системой DV-Geo, не все о ней знают. Помимо данного курса есть подробная документация, и все необходимые вопросы можно освоить самостоятельно (или с помощью коллег).

Важно то, что мы пока ничего не знаем о том, как строится модель. Методику построения геологической модели (на примере модели «CGE») мы рассмотрим в Части II.

63

Часть II. Построение геологической модели «CGE»

Мы начинаем вторую часть учебного курса по DV-Geo, в которой мы воспроизведем процесс построения только что рассмотренной модели «CGE» отталкиваясь от ее исходных данных. Сразу же возникает важный вопрос - какие данные мы считаем исходными? Для DV-Geo это, как правило, координаты устьев скважин на плане XY, данные инклинометрии, альтитуды, данные ГИС, контуры разломов. В отношении исходных данных именно модели «CGE» необходимо обратить внимание на два обстоятельства. Первое - мы видели, что в рассматриваемой модели используется вертикальная аппроксимация скважин. Поэтому вместо данных инклинометрии у нас будут координаты пластопересечений скважин и целевого горизонта на плане XY. Названное обстоятельство не является существенным, поскольку ничего в процессе построения модели не меняет. Второе же обстоятельство, напротив, является очень существенным. Поскольку работа по созданию учебной модели все же не должна требовать от пользователя столь больших затрат времени и усилий, как это бывает при исполнении реального проекта, набор действительно первичных данных дополнен готовыми результатами корреляции скважин, колонками РГИ и отметками ВНК. То есть, два важных этапа построения геологической модели в DV-Geo (корреляцию скважин и интерпретацию ГИС) можно при желании пропустить. Все названные исходные данные находятся в каталоге “Geodat\testdata”.

64

1. Создание нового проекта. Загрузка исходных данных

1.1. Создание нового проекта

Запустим систему двойным щелчком левой клавишей мыши по ее ярлыку на рабочем столе. Откроется главное окно DV-Geo, которое уже было показано на фиг. 6. Вместо команды «Файл\Загрузить» выберем в меню главного окна команду «Файл\Сохранить как» и введем имя «MyCge». Обратите внимание, сохранение должно производиться в каталоге «Geodat». В строке «Тип файла» при этом должно быть указано «Проект (создание нового проекта)» (фиг. 52).

Фиг. 52. Сохранение нового проекта под именем «MyCge»

После нажатия на кнопку «Сохранить» происходит следующее. Первое - в каталоге «Geodat» появляется папка «MyCge» со стандартным набором вложенных каталогов для хранения файлов модели. Пока все эти вложенные каталоги пусты. Второе – в каталоге «Geodat» появляется файл «MyCge.dvg», который мы должны будем загружать в начале каждого сеанса работы с этой моделью. Третье – имя «MyCge.dvg» появляется в заголовке главного окна DVGeo.

Окно «Список объектов» нового проекта открывается командой главного меню «Данные\Загруженные данные», и далее через меню пункта «Current» в окне «Загруженные данные» (фиг. 53). Можно убедиться, что в начальный момент все группы данного окна (за исключением групп «Tools», «Patterns» и «Environment») пусты.

65

Фиг. 53. Открытие окна «Список объектов» проекта «MyCge»

1.2. Загрузка траекторий скважин и альтитуд

Первыми данными, загружаемыми в модель, являются траектории скважин. Для загрузки скважин используется процедура «Загрузить скважины», написанная на языке Th и расположенная в группе «Tools» (фиг. 54). Подробное описание загрузки скважин дано в томе 1 «Создание проекта». Еще более подробное описание, включающее построковый разбор текста процедуры на языке Th, представлено в Дополнительном томе 1 «Организация данных и процедур».

Фиг. 54. Процедура загрузки траекторий скважин

66

Поэтому мы можем быть краткими. Процедура «Загрузить скважины» запускается через ее собственное меню (командой «Выполнить»). Сразу после запуска процедура запрашивает файл с координатами устьев (или пластопересечений) скважин на плане XY. Предполагается, что этот файл имеет имя «coord.txt», но возможно и другое имя. Файл с координатами должен иметь текстовый формат, показанный на фиг. 55.

Фиг. 55. Формат файла с координатами устьев скважин

Файл «coord.txt» с координатами скважин для модели «CGE» находится в каталоге «Geodat\Testdata». Получив имя этого файла, процедура загрузки открывает диалоговое окно, показанное на фиг. 56.

Фиг. 56. Диалоговое окно процедуры загрузки скважин

Поясним (тоже коротко) значение основных полей. В слотах «Group» и «Cluster» указываются имена группы и подгруппы в окне «Список объектов», в которые будут помещены

67

загруженные скважины. Имя группы должно быть «SKV», имя подгруппы может быть произвольным (8 символов, допускаются буквы, цифры, и символы подчеркивания). В поле «Status» указывается класс скважин. Здесь, наряду с универсальным значением «SKV», возможны более конкретные значения «OIL», «WAT», «OIL_WAT» и некоторые другие. Поля «Mult X» и «Mult Y» позволяют умножить значения координат в файле на заданные коэффициенты. Значение «Load» в поле «Space Bounds» означает, что границы пространства XY будут установлены по минимальным и максимальным координатам загружаемых скважин, с небольшим запасом, определяемым значениями в полях «Zone X» и «Zone Y».

Важным является значение в поле «No LST files». Процедура устроена так, что, загрузив координаты XY очередной скважины, она начинает искать в том же каталоге стандартный LST-файл (с расширением «*.lst»), имя которого совпадает с именем скважины. Если искомый файл не находится, процедура может или пропустить эту скважину, или ввести ее как вертикальную, в зависимости от значения в поле «No LST files».

В модели «CGE» все скважины в пределах целевого горизонта считаются вертикальными. При этом в файле «coord.txt» стоят координаты не устьев скважин, а пластопересечений. Это позволяет нам обойтись без LST-файлов. Соответственно, при загрузке скважин в наш проект в поле «No LST files» должно стоять значение «Считать скважину вертикальной».

Выбранное значение в последнем поле «Load» определяет действия процедуры в случае, когда в модели уже есть скважина, одноименная c загружаемой. Опция «Все» означает, что старую скважину следует заменить на новую. Опция «Только новые» означает, что одноименную скважину повторно грузить не следует.

В случае успешной загрузки скважины появляются в группе «SKV\Well» (или в иной указанной группе) окна «Список объектов». Полностью убедиться в успешной загрузке скважин мы можем после того, как бросим их на план XY (фиг. 57). Обратите внимание, границы нашего пространства XY ровно на 3000 м отстоят от крайних скважин.

Фиг. 57. Результат загрузки скважин в проект «MyCge»

68

Загрузка альтитуд скважин выполняется при помощи Th-процедуры «Tools\Converters\LST\Загрузить поправку». После запуска процедуры мы должны указать имя файла с поправками. Для модели «CGE» файл с альтитудами «popr.alt» также находится в каталоге ««Geodat\Testdata». Формат файла с альтитудами показан на фиг. 58.

Фиг. 58. Формат файла с альтитудами

Значения альтитуд скважин можно увидеть в окне «Профиль» при отображении скважин в режиме «Простой», если для фреймов скважин выбрана опция «Показывать альтитуду» (фиг. 59). Альтитуды отображаются под именами скважин.

Фиг. 59. Результат загрузки альтитуд скважин в окне «Профиль»

Мы обошлись при загрузке скважин без LST-файлов. Тем не менее, отметим, что их формат подробно описан и в Томе 1, и в Дополнительном томе 1.

69

После загрузки скважин необходимо выполнить сохранение модели. При этом в каталоге «MyCge\default\LST\Well» появится бинарный файл внутреннего формата DV-Geo с расширением «wells.vfs». Этот файл содержит таблицы инклинометрии скважин. Координаты скважин на плане XY и соответствующие альтитуды содержит соседний текстовый файл «xyz.txt». Список с именами скважин появятся также в файле «MyCge.dvg».

1.3. Загрузка данных ГИС

На фиг. 57 можно видеть, что в окне “Список объектов” в группах с именами скважин присутствуют только траектории скважин и таблицы “Конструкция”. Для ввода в модель данных ГИС необходимо вызвать Th-процедуру «Tools\Converters\LST\Загрузить объекты Gis». Подробное описание работы с процедурой дано в томе 1 «Создание проекта». Еще более подробное описание, включающее разбор текста процедуры на языке Th, представлено в дополнительном томе 1 «Организация данных и процедур». Поэтому ограничимся кратким описанием.

Запуск процедуры «Загрузить объекты Gis» производится командой «Выполнить» из ее меню. После запуска открывается диалоговое окно (фиг. 60).

Фиг. 60. Первое диалоговое окно процедуры загрузки данных ГИС

Значение в первом поле «Из всех» означает, что перед загрузкой процедура просмотрит заголовки всех загружаемых LAS-файлов и предложит пользователю выбрать нужные методы ГИС из полного списка методов, представленных в этих файлах. Альтернативное значение «Набор по умолчанию» исключает предварительный просмотр заголовков. В этом случае процедура предложит пользователю выбирать методы из фиксированного стандартного списка «PS, BK, NKT, IK, GK». Во втором поле указывается расширение LAS-файлов. В третьем поле указывается кодировка русского алфавита, необходимая для корректного считывания из LAS-файлов названий методов ГИС на русском языке.

После нажатия кнопки «ОК» процедура открывает стандартное диалоговое окно «Открыть файл», в котором надо указать загружаемый LAS-файл. Одновременно в этом диалоге фиксируется вмещающий каталог. В нашем примере необходимо указать файл «Geodat\Teatdata\Las\157.las». После ввода имени файла открывается третье диалоговое окно процедуры, представленное на фиг. 61.

70

Фиг. 61. Третье диалоговое окно процедуры загрузки данных ГИС

В данном окне мы определяем, для какой группы скважин производится загрузка данных ГИС, в какую группу будут помещены эти методы ГИС на каждой скважине и какие именно методы надо загрузить (из стандартного списка или из полного). Также мы указываем режим добавления или замещения кривой в случае, если одноименная кривая на скважине уже имеется. В режиме «Добавить» при совпадении имен новая кривая будет получать дополнительный порядковый номер в скобках, например GK, GK(1), GK(2) и т.д.

Поле «Имена скважин» позволяет указать, какими могут быть имена LAS-файлов до расширения. Эти имена должны или повторять имена соответствующих скважин, или содержать имена скважин как подстроку. Третий вариант позволяет читать данные каротажа только из одного явно указанного файла. Поля «Верхняя/Нижняя граница» позволяют задать интервал загрузки данных ГИС.

Ход загрузки кривых отображается в окне «Выполнение Th». По окончании процедуры открывается окно программы «Блокнот» со списком не найденных LAS-файлов для скважин указанной группы. Загруженные каротажные кривые начинают отображаться в окне «Список объектов» как элементы группы «Gis» (если было указано именно это имя) своих скважин. Окончательно убедиться в успешной загрузке кривых можно после их отображения в окне «Профиль» (фиг. 62).

71

Фиг. 62. Отображение загруженных каротажных кривых в окне «Профиль»

После загрузки кривых необходимо выполнить сохранение модели. При этом в каталоге «MyCge\default\LAS\Well» появятся бинарные файлы внутреннего формата DV-Geo с расширением «* .lfs». Соответствующие ссылки появятся и в файле «MyCge.dvg».

1.4. Загрузка данных РГИ

Для ввода в модель данных РГИ необходимо вызвать Th-процедуру «Tools\Converters\LST\Загрузить объекты Rgi». Описание работы с процедурой дано в томе 1 «Создание проекта». Поэтому ограничимся кратким изложением.

Запуск процедуры «Загрузить объекты Gis» производится командой «Выполнить» из ее меню. При этом открывается диалоговое окно (фиг. 63).

Фиг. 63.

Первые два поля аналогичны полям на фиг. 60. Числа в третьем и четвертом поле позволяют настраиваться на переменные форматы таблиц Rgi. Первое число указывает количество блоков заголовка (отмеченных знаком «~»), которые предшествуют списку параметров. Второе число указывает количество колонок, которые предшествуют первой колонке с

72

параметром. Маркеры секции скважин и секции кривых показывают символы (два варианта), после которых следуют загружаемые данные. Для загрузки данных из приложенных файлов описанные поля должны быть заполнены по образцу (причем «С» должны быть буквами русского алфавита). После нажатия кнопки «ОК» открывается стандартный диалог «Открыть файл», в котором необходимо указать каталог и первый файл с данными Rgi. В нашем случае это должно быть «Geodat\Teatdata\Rgi\157.rgi». Затем открывается третье диалоговое окно, показанное на фиг. 64.

Фиг.64.

Значения в полях этого окна совпадают по смыслу с теми, которые запрашивала процедура загрузки данных ГИС (см. фиг. 61).

Ход загрузки кривых Rgi отображается в окне «Выполнение Th». В случае успешной загрузки кривые Rgi начинают отображаются в окне «Список объектов» как элементы группы «Rgi» (если было указано именно это имя) своих скважин. Окончательно убедиться в успешной загрузке кривых можно после их отображения в окне «Профиль».

Поскольку процедура загрузки данных Rgi написана на языке Th, ее можно настроить на любую модификацию формата входных данных. После загрузки кривых необходимо выполнить сохранение модели. При этом в каталоге «MyCge\default\RGI\Well» появятся бинарные файлы внутреннего формата DV-Geo с расширением «* .rfs». Соответствующие ссылки появятся и в файле «MyCge.dvg».

В отличие от данных ГИС, колонки РГИ являются не внешними данными, а результатом интерпретации каротажных кривых в системе DV-Geo. В настоящем учебном курсе мы не касаемся методик интерпретации данных ГИС. Используемые при этом инструменты описаны в томе 6 «Обработка скважинных данных».

1.5. Загрузка корреляции скважин

Готовые результаты корреляции скважин нашего примера загружаются в модель командой «Средства\Импорт\Загрузить» из меню группы «Correlation» (фиг. 65). Обратите внимание, из состава меню видно, что в DV-Geo реализованы процедуры загрузки данных корреляции из трех внешних форматов. Напомним также, что процедуры, вызываемые через пункт меню

73

«Средства», написаны на языке Th, и их легко можно модифицировать под любой внешний формат.

Фиг. 65. Вызов процедуры загрузки данных корреляции

Сразу после запуска процедуры открывается диалоговое окно «Открыть файл», в котором необходимо указать имя файла «Geodat\Teatdata\Correlation.crr». Далее открывается второе диалоговое окно, показанное на фиг. 66. В поле «Глубина» необходимо указать кабельную длину «L», поскольку числовые значения в колонках файла «Correlation.crr» показывают именно кабельную длину.

Фиг. 66. Диалог загрузки данных корреляции

В случае успешной загрузки корреляционные отметки появляются в группе «Correlation» Текстовый формат файла «Correlation.crr» показан на фиг. 67.

Фиг. 67. Формат файла с данными корреляции.

74

Загрузка скважинных отметок уровня ВНК из файла «OWC7_1.TXT» производится аналогично.

1.6. Загрузка объектов «Area»

Загрузим линию разлома из файла «Geodat\Teatdata\Fault.txt». Сначала создадим для нее группу «Разлом» внутри группы «Area». Это делается через меню группы «Area». Затем через меню группы «Разлом» вызовем процедуру «Средства\Импорт\Добавить линию (DV-Geo)» (фиг. 68).

Фиг. 68. Загрузка линии разлома

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

Фиг. 69. Диалоговое окно загрузки объекта «Area»

В результате исполнения процедуры линия разлома появляется в группе «Разлом», откуда ее можно бросить на план XY. Формат внешнего файла с координатами точек разлома показан на фиг. 70.

75

Фиг. 70. Формат внешнего файла с координатами точек разлома

1.7. Загрузка данных перфорации

Данные перфорации скважин загружаются в модель Th-процедурой «Tools\Converters\LST\Загрузить перфорацию». Описание работы с процедурой дано в томе 1 «Создание проекта». Поэтому ограничимся кратким пояснением.

Запуск процедуры производится командой «Выполнить» из ее меню. После этого открывается диалоговое окно (фиг. 71), в котором необходимо указать имя, под которым данные перфорации будут помещены в группу «Perforation» каждой скважины. Имя группы «Perforation» менять не следует. Далее потребуется указать файл с данными перфорации (в нашем случае это файл «Geodat\Teatdata\Perf.prf»). После этого загрузка будет произведена. Для отображения данных перфорации той или иной скважины необходимо перенести элемент «Perf» из ее списка объектов на ее же отображение в окне «Профиль» (фиг. 72). Примечание: элемент «Perf» появляется в списке объектов всех скважин, в том числе и тех, для которых никаких данных перфорации нет.

Фиг. 71. Диалоговое окно загрузки данных о перфорации

76

Фиг. 72. Данные перфорации скважины «327» в окне «Профиль», в окне «Список объектов», и в исходном текстовом файле.

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

77

2. Построение структурной модели

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

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

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

В DV-Geo мы можем выполнить и ручную, и автоматическую корреляцию скважин. Кроме того, мы можем (и мы это уже рассматривали) загрузить в DV-Geo результаты корреляции из внешних файлов. Причем одно здесь не исключает другое. Когда мы делаем корреляцию сами, мы всегда начинаем с того, что загружаем в модель все варианты корреляции, выполненной ранее.

2.1. Ручная корреляция скважин по данным ГИС

Мы уже знаем, что основным инструментом при выполнении ручной корреляции скважин является окно «Профиль», ассоциированное с прямоугольным селектором скважин на плане XY. Создание названного селектора и открытие окна «Профиль» мы рассматривали в Части I, в разделе 5. Мы также знаем, как ставятся отметки корреляционного репера. Напомним последовательность действий при ручной корреляции:

1.Через меню группы «Correlation» создается новый репер, имя которого должно быть оригинальным. Длина имени не должна превышать 8 символов, допускаются буквы, цифры и знак подчеркивания (фиг. 73).

78

Фиг. 73. Создание нового репера

2.Вновь созданный репер переносится из списка объектов на рабочую область окна «Профиль». При этом его имя появляется в «выпадающем списке» на панели окна «Профиль» (фиг. 74). Отметок на скважинах (если это вновь созданный, а не корректируемый репер) пока еще нет.

3.Режим установки (корректировки) репера устанавливается нажатием кнопки «Двигать метку» (фиг. 74).

4.Отметка репера устанавливается нажатием левой кнопки (и последующим небольшим смещением) мыши на фрейме скважины. Вариант отображения фреймов («Полный», «Простой» и т.д.) при этом может быть любым.

5.Перейти к установке другого репера (при условии, что он тоже был брошен в окно «Профиль») можно или выбрав его имя в выпадающем списке, или захватив одну из его отметок нажатием левой кнопки мыши. В последнем случае переключение на другой устанавливаемый репер производится автоматически, и за этим необходимо следить. Чтобы исключить случайные смещения отметок ранее установленного репера, желательно его заблокировать (через меню репера в окне «Профиль»).

6.В процессе корреляции возможно отображение установленных реперов на плане XY. Описание данной опции приведено в Части I, в разделе 5.3.

Установка отметок корреляционного репера, как правило, производится в условиях выравнивания скважин по вышележащему реперу, установленному или загруженному ранее. Названное выравнивание, напомним, производится через меню репера в окне «Профиль». На фиг. 74 выравнивание произведено по реперу «top71c». Если репер, по которому производится выравнивание, установлен не на всех скважинах, визуальный эффект выравнивания пропадает. Поэтому рекомендуется сочетать выравнивание скважин с фильтрацией скважин на плане XY по реперу, задающему выравнивание. Как уже было сказано, возможно также выравнивание скважин по двум ранее установленным реперам.

79

Фиг. 74. Установка отметок корреляционного репера

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

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

Более подробно инструментарий ручной корреляции описан в томе 3-1 «Корреляция разрезов скважин».

На исследуемой площади была выполнена ручная корреляция скважин по данным ГИС. В результате на скважинах были выделены два стратиграфических горизонта: «71с» (кровлей является репер «top71c», подошвой - «bot71c») и «72с» (соответственно реперы «top72c» и «bot72c»). В рамках описываемого проекта моделируется только стратиграфический горизонт

«71с».

80

2.2. Автоматическая корреляция скважин по данным ГИС

Поскольку автоматическая корреляция скважин при построении модели «CGE» не использовалась, опишем ее максимально кратко. В основе автоматической корреляции скважин лежит понятие «квазистратиграфическая поверхность». Что это такое, поясняется на фиг. 75.

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

Вверхней части фиг. 75 изображен условный геологический разрез, стратиграфические горизонты которого частично размыты. Из представленного рисунка следует, что задача автоматического прослеживания на скважинах конкретной стратиграфической поверхности невыполнима, поскольку объективные основания для такого прослеживания часто (иногда на большинстве скважин) отсутствуют.

Всредней части фиг. 75 оставлены в цвете только те горизонты предыдущего разреза, на которых корреляция между соседними скважинами имеет объективные основания. Из этой части иллюстрации следует, что если рассматривать не одну стратиграфическую поверхность, а некоторый протяженный стратиграфический интервал (мощностью 50-60 м), то в его пределах объективные основания для корреляции скважин (вверху, в середине, или внизу интервала) существуют всегда. То есть, мы можем рассчитать некоторую «квазистратиграфическую» поверхность, ассоциированную с центром этого интервала, и идущую параллельно рельефу наиболее выразительного (в пределах этого интервала) ненарушенного горизонта (фиг. 75 внизу). Понятно, почему мы называем эту поверхность квазистратиграфической - она близка к настоящей стратиграфической поверхности там, где последняя существует, и разумно интерполирует ее в областях размывов.

Уточним следующее. Автоматическая корреляция скважин дает нам не саму квазистратиграфическую поверхность, а только ее отметки на скважинах (автоматический

81

репер). Многократное выполнение процедуры автоматической корреляции дает нам набор таких реперов, смещенных друг относительно друга, допустим, на полуширину расчетного интервала, то есть на 25-30 м. Можно считать эти репера и через 10 м, но при ширине расчетного интервала в 50-60 м большого смысла в этом не будет. Интерполяция отметок этих реперов дает нам стопку квазистратиграфических поверхностей, то есть каркас модели. Но напомним, что на этапе интерполяции мы должны будем учесть рельеф сейсмических отражающих поверхностей и поверхности тектонических нарушений (если они у нас будут).

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

1.Обнаруживаются все ошибки традиционной ручной корреляции, которых больше, чем принято считать.

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

Вместе с тем автоматическая корреляция требует стандартизации данных ГИС. Различия каротажных кривых должны быть обусловлены только различиями геологической среды, но не условиями наблюдения, типом аппаратуры и т.д. При использовании автоматической корреляции основную часть времени занимает этап предварительного редактирования данных. Опыт показывает, что наиболее эффективной автоматическая корреляция оказывается при выполнении работ на площадях с числом скважин сто и более. Подробное описание алгоритма автоматической корреляции, методики ее применения и примеров использования дается в томе 3-2 «Корреляция разрезов скважин с использованием расширенных средств поддержки».

2.3. Создание региона

В проекте «CGE» каркас модели включает только две поверхности - кровлю и подошву горизонта «71с». Названные поверхности получены посредством интерполяции соответствующих отметок на скважинах. Первым шагом в расчете этих поверхностей (если не считать этап корреляции данных ГИС) является создание региона – сетки на плане XY. Мы достаточно подробно рассматривали регион и порядок его создания при знакомстве с моделью в Части I. Поэтому сейчас мы знаем, что необходимо сделать. Первое - создать регион через меню группы «Region» окна «Список объектов», второе – кинуть его на план XY и откорректировать его границы мышью, третье - установить (через меню региона на том же плане) размер ячеек его сетки равный 50х50 м.

Что произойдет, если координаты региона в нашей модели будут несколько отличаться от координат региона в модели «CGE»? Это первый пункт, в котором мы можем нечаянно отклониться от нашего образца. Ничего страшного, небольшие отличия не страшны. Важно только, чтобы наш регион охватил все скважины и чтобы он охватил будущий внешний контур ВНК. При выполнении двух названных условий объемы запасов в нашей новой модели и в нашем образце совпадут. О кратности размеров региона размеру ячеек можно не беспокоиться, поскольку число ячеек сетки по каждой стороне региона округляется до целого в большую сторону (и фактические размеры региона чуть-чуть увеличиваются).

82

2.4. Интерполяция структурных поверхностей

После того, как регион «1» создан, мы можем рассчитать на нем (то есть используя его сетку) нужные нам поверхности. Расчет поверхности выполняется в два шага. Первый шаг – в меню региона выбирается пункт «Добавить поверхность». Дальше в диалоговом окне необходимо определить имя новой поверхности («1») и выбрать из списка ее скважинный репер (фиг. 76). Значения Nx и Ny по умолчанию соответствуют сетке региона, менять их не следует.

Фиг. 76. Создание новой поверхности

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

Второй шаг создания поверхности – в ее меню выбирается пункт «Вычислить». В результате открывается диалоговое окно со многими вкладками (фиг. 77). В простейшем случае интерполяции поверхности по одному только скважинному реперу нам нужна в этом окне только одна вкладка – «Локальная интерполяция» (она и показана). На ней мы указываем следующее:

•в поле «Операция» - должна быть произведена именно интерполяция поверхности, а не локальное уточнение поверхности по новым скважинам;

•опция «Сечение» (в соседнем поле) означает, что координаты X, Y, Z отметок репера должны быть рассчитаны по данным инклинометрии. Альтернативное значение «План» указывается в случае, когда интерполируемая координата Z берется непосредственно из таблицы значений репера (или объекта), а недостающие X, Y берутся из файла «coord.txt»;

•в поле «Репер» мы выбираем имя интерполируемого репера. В следующем поле «Объект» мы ничего не выбираем. Здесь мы могли бы выбрать имя таблицы «скважина - параметр», но два последних поля альтернативны друг-другу;

•в двух первых полях параметра «Вес» указывается показатель степени весовой функции 1/Rn, отдельно для направлений X и Y. В третьем поле на этой горизонтали вводится предельное число ближайших скважин, влияющих на интерполяцию в текущей точке. В этих трех полях лучше ничего не менять;

•в поле «Радиус» указываются предельные радиусы интерполяции для минимальной и максимальной плотности скважин. Фактическое значение радиуса определяется на основе локальной оценки плотности скважин. Два значения в следующей строке позволяют

83

придать зоне влияния эллиптическую форму. Первое из них задает отношение осей эллипса, второе - угол наклона главной оси. Мы анизотропию не задаем, и значения радиусов «по умолчанию» приемлемы;

•сглаживание не производится;

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

Остальные параметры могут также иметь значения «по умолчанию». Полное описание параметров интерполяции поверхности дано в томе 4 «Построение трехмерной структурной модели».

Фиг. 77. Диалоговое окно расчета поверхности

После заполнения вкладки «Локальная интерполяция» следует нажать на ней кнопку «Вычислить». При этом открывается окно процедуры, в котором будет отображаться ход вычислений. Во время расчета поверхности диалог с системой DV-Geo не блокируется, то есть пользователь может продолжать работу с проектом. В случае успешного завершения расчетной процедуры имя поверхности в списке объектов начинает отображаться обычным цветом.

Если мы планируем использовать рассчитанную поверхность не только в текущем, но и в следующих сеансах работы с проектом, необходимо сохранить модель. При этом в каталоге

84

«MyCge\default\CORRL\surf» появится файл поверхности во внутреннем (текстовом) формате DV-Geo с расширением «*.crs». Соответствующая ссылка появится и в файле

«MyCge.dvg».

Поверхность «2» (по реперу «bot71c» ) рассчитывается аналогично. Но лучше внести одно небольшое изменение. А именно, выбрать пункт «Добавить поверхность» не в меню региона, а в меню поверхности «1». Разница будет в том, что все параметры новой поверхности будут «по умолчанию» копировать параметры поверхности «1».

Для оценки качества рассчитанных поверхностей необходимо визуализировать их в окне 3D, на плане XY и в окне «Профиль». Необходимо убедиться, что поверхности не имеют «выбросов», разрывов и что они точно проходят через отметки на скважинах. Обе рассчитанные поверхности показаны на фиг. 78.

Фиг. 78. Качество структурных поверхностей «1» и «2» контролируется в окнах «План слайса» и «Профиль»

На этом в проекте «CGE» построение структурной модели завершается. Но мы еще один раз используем инструмент интерполяции поверхностей для расчета поверхности ВНК.

2.5. Интерполяция поверхности ВНК

Поверхность водонефтяного контакта рассчитывается по скважинному реперу «OWC7_1» с теми же параметрами на вкладке «Локальная интерполяция», что и структурные поверхности. Полученная карта поверхности ВНК показана на фиг. 79. Можно видеть, что скважина 158Р

85

задает аномально низкий уровень ВНК. Соответствующую скважинную отметку желательно перепроверить.

Фиг. 79. Поверхность ВНК, полученная посредством интерполяции отметок ВНК на скважинах

86

3. Литологическое моделирование

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

Куб литологии должен (первое) воспроизводить литологические колонки скважин и (второе) соответствовать принципиальной модели геологической среды. Принципиальной моделью называют карту эффективной мощности коллектора, которую геолог строит на основе скважинных данных и своих представлений о генезисе исследуемого объекта. Для получения именно такого куба предлагается придерживаться следующей последовательности действий:

1.Посредством интерполяции кривых РГИ параметра альфа-ПС на скважинах рассчитывается куб «aps» (для ограниченного применения).

2.Из куба «aps» при помощи калькулятора рассчитывается индексный куб литологии. Если значение aps в ячейке куба больше 0.35, то в ячейку куба литологии заносится 1 (коллектор), в противном случае заносится 0 (неколлектор). Пороговое значение aps может корректироваться.

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

4.В силу более грубой вертикальной дискретизации трехмерной модели в сравнении со скважинными данными РГИ куб литологии завышает эффективные мощности в точках скважин. Для устранения названной погрешности используется куб NTG. В кубе NTG содержатся интерполированные скважинные значения параметра «доля коллектора в ячейке куба».

Далее вышеперечисленные этапы литологического моделирования описываются подробно.

3.1. Расчет куба APS

Выберем в меню группы «Pack» пункт «Добавить куб» (фиг. 80). В появившемся диалоговом окне зададим имя куба («aps») и укажем имя региона. Должен быть указан именно тот регион, на котором рассчитаны структурные поверхности модели. После этого имя нового куба появится в группе «Pack».

Фиг. 80. Добавление куба «aps» в группу «Pack»

87

Если мы уже имеем кубы в группе «Pack», то пункт «Добавить куб» можно вызвать из меню одного из них. В этом случае все исходные параметры нового куба будут воспроизводить параметры куба, которому принадлежало меню.

Из меню добавленного куба «aps» (пока его имя отображается в группе «Pack» бледным цветом) выберем пункт «Вычислить». После этого открывается диалоговое окно расчета куба со многими вкладками. Две из них, «Исходные данные» и «Локальная интерполяция» (другие вкладки нам не потребуются) показаны на фиг. 81.

Фиг. 81. Диалоговые вкладки с параметрами расчета куба «aps»

Для расчета куба «aps» значения параметров должны быть установлены так, как показано на иллюстрации. Смысл большинства параметров на вкладке «Исходные данные» достаточно очевиден. Мы должны отметить поверхности структурного каркаса, указать исходные данные Rgi.Апс, которые будут интерполироваться, и уточнить Nz - количество слоев ячеек стратиграфической сетки по Z. Значения Nx и Ny менять не следует, они воспроизводят значения сетки региона.

Далее, в строке «Разделять» первой вкладки следует указать «Rgi.Литология» и в соседнем поле поставить значение «1». Смысл этого в следующем. Поскольку дискретизация скважинных кривых РГИ, как правило, меньше дискретизации куба по вертикали, при заполнении ячеек куба, лежащих на траекториях скважин, происходит осреднение по нескольким значениям Rgi.Апс. Но из раздела «Введение» (пункт 1.1) мы знаем, что смешивать при осреднении значения коллектора и неколлектора нежелательно. Параметры в строке «Разделять» определяют, что если, судя по кривой Rgi.Литология, в ячейку куба попадают коллектор и неколлектор, осреднение Rgi.Апс производится только для значений, соответствующих коллектору.

88

Но такое «правильное» осреднение имеет побочное следствие - эффективная мощность коллектора в кубе, по сравнению с данными скважин, оказывается завышенной. Действительно, при попадании в ячейку хотя бы одного значения РГИ «коллектор», вся ячейка получает значение Апс коллектора. При этом ошибка тем больше, чем больше расчлененность разреза. Для устранения названной ошибки приходится использовать дополнительный куб - куб NTG. Расчет этого куба будет описан ниже.

На вкладке «Локальная интерполяция» (необходимо применять именно этот вид интерполяции) мы указываем следующее:

•должна быть произведена именно интерполяция (а не иные предусмотренные на вкладке операции);

•интерполяция должна производиться согласно структуре, то есть по слоям стратиграфической сетки (а не горизонтально, как, допустим, насыщенность);

•исходными для интерполяции являются значения в ячейках, полученные посредством взвешенного осреднения данных РГИ в скважинах (с разделением по литологии);

•два предельных радиуса интерполяции, один для области с максимальной плотностью скважин, другой – для области с минимальной плотностью скважин. В точке интерполяции рассчитывается локальная плотность скважин, после чего пропорционально определяется радиус;

•сглаживание не производится;

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

Остальные параметры могут иметь значения «по умолчанию». Полное описание параметров интерполяции кубов дано в томе 5 «Построение многопараметровой модели». После заполнения вкладок «Исходные данные» и «Локальная интерполяция» следует нажать кнопку «Вычислить» на второй из них. При этом открывается окно процедуры, в котором отображаются послойные (по каждому слою ячеек стратиграфической сетки) вычисления куба. Информация, выводимая в окно процедуры, позволяет прогнозировать время расчета, которое может быть значительным. Во время расчета куба диалог с системой DV-Geo не блокируется, то есть пользователь может продолжать работу с проектом. Рассчитанный куб «aps» вместе с исходными кривыми РГИ показан на фиг. 82.

Не противоречит ли расчет куба «aps» нашим предыдущим утверждениям о том, что скважинные данные нельзя интерполировать? Нет, не противоречит. Во-первых, мы не собираемся использовать куб «aps» для прогноза пористости. Он нам нужен только для расчета куба литологии. Во-вторых, полученный куб литологии мы еще будем подгонять и под принципиальную модель, и под скважинные данные.

89

Фиг. 82. Исходные кривые РГИ «Aпс» и рассчитанный по ним куб «aps» в окне «Профиль»

3.2. Расчет куба литологии

Расчет индексного куба литологии «lito» производится из куба «aps» при помощи калькулятора. Значение «aps» большее 0.35 интерпретируется как «коллектор» (и в соответствующую ячейку куба литологии заносится 1), меньшее 0.35 – как неколлектор (и в ячейку куба литологии заносится 0). Пороговое значение 0.35 может корректироваться. Вид окна калькулятора с условиями расчета и полученный куб литологии показаны на фиг. 83, 84.

Фиг. 83. Окно калькулятора для расчета куба литологии по кубу «aps»

90

Фиг. 84. Куб «lito» в окне «Профиль». Обратите внимание - эффективные мощности в кубе получаются завышенными в сравнении с данными «Aпс» по скважинам

Для оценки куба литологии необходимо рассчитать соответствующую ему карту эффективной мощности коллектора. С этой целью в меню куба «lito» должен быть активирован пункт «Средства\Calculate Map\Эффективная мощность» (фиг.85). Появляющееся после этого диалоговое окно, установленные в нем параметры и результат показаны на фиг. 86.

Фиг. 85. Расчет карты эффективной мощности по кубу «lito»

91

Фиг. 86. Расчет карты эффективной мощности по кубу литологии без сглаживания

Карта эффективной мощности показывает, что куб «lito» не является удовлетворительным по двум причинам. Первая причина – изменение мощности коллекторских тел в пространстве не соответствует принципиальной модели геолога. Вторая причина - мощности песчаных тел в точках скважин несколько превышают мощности, рассчитанные по самим скважинам. Для устранения названных недостатков выполняют (первое) коррекцию куба «lito» по принципиальной модели и (второе) расчет дополнительного куба «NTG».

3.3. Коррекция куба литологии по принципиальной модели

Выполним следующий численный эксперимент. Рассчитаем по кубу литологии вторую карту эффективной мощности «HefPr» с тем отличием, что применим сильное сглаживание (Smooth) и сильное подтягивание (Drag).

Назовем полученный результат (фиг. 87) картой принципиальной модели, за неимением действительно таковой. Это возможно, поскольку карта «HefPr», благодаря сглаживанию, выглядит достаточно привлекательно. В то же время, эта карта (снова благодаря сглаживанию) уже не соответствует кубу «lito». Покажем, как можно использовать карту принципиальной модели для коррекции куба «lito».

92

Фиг. 87. Расчет карты эффективной мощности «HefPr» по кубу литологии при сильном сглаживании и сильном подтягивании

Инструментами решения названной задачи являются шаблоны из группы «Environment\Map3D» окна «Список объектов». Создадим в упомянутой группе (через ее меню) новый шаблон, который назовем «Principal» (фиг. 88).

Фиг. 88. Шаблон «Principal» для коррекции куба литологии

Откроем шаблон «Principal» (уже через его собственное меню) для редактирования. Появится окно «Map3D.Principal» без текста, разделенное на две неравные части. Щелкнем два раза левой кнопкой мыши в левой части окна. Откроется диалоговое окно для ввода имени куба. Введем имя «LitoPrinc», как показано на фиг. 89.

93

Фиг. 89. Окно «Map3D.Principal». Ввод имени куба

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

Фиг. 90. Заполнение полей шаблона для коррекции куба «lito» по карте принципиальной модели «HefPr»

Нажатие кнопки «Литология» запускает расчет куба «LitoPr». Карта эффективной мощности куба «LitoPr» будет без сглаживания и подтягивания воспроизводить основные черты карты принципиальной модели. Сравнение исходного и модифицированного куба литологии в окне «Профиль» показано на фиг. 91.

94

Фиг. 91. Сравнение кубов «lito» (вверху) и «LitoPr» (внизу). Карта эффективной мощности куба «LitoPr» приближается к карте эффективной мощности принципиальной модели

Вместе с тем, карта эффективной мощности куба «LitoPr» по-прежнему не может точно воспроизвести эффективные мощности по скважинным данным, в силу того, что размер ячейки куба по Z существенно больше интервала дискретизации кривых РГИ. Для решения этой проблемы используется куб NTG (от английского «net to gross»).

3.4.Расчет куба NTG

Вячейках куба NTG на траекториях скважин стоят значения, равные отношению мощности коллектора в пределах ячейки по скважине к мощности ячейки. Весь куб NTG рассчитывается посредством интерполяции значений на скважинах. При всех последующих вычислениях (карт эффективных мощностей, карт пористости, объема запасов) высота ячейки коллектора в кубе «LitoPr» будет умножаться на значение из соответствующей ячейки куба NTG.

Параметры расчета куба NTG показаны на фиг. 92. Параметры на первой вкладке те же, что и при расчете куба «aps». На второй вкладке тип осреднения данных в ячейке установлен как «Отношение». Кроме того, используется логический фильтр по кубу литологии. Вид рассчитанного куба NTG в окне «Профиль» показан на фиг. 93. На этом литологическое моделирование завершается.

95

Фиг. 92. Диалоговые вкладки с параметрами расчета куба «NTG»

Фиг. 93. Куб NTG в окне «Профиль»

96

4. Петрофизическое моделирование

4.1. Расчет куба пористости

Итак, куб литологии и куб NTG у нас есть. Расчет куба пористости производится посредством интерполяции данных РГИ о пористости, но только в пределах коллектора. При этом наличие в пространстве зон неколлектора игнорируется. Этот подход уязвим для критики, но другого способа расчета пористости в категориальной среде у нас нет.

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

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

Фиг. 94. Диалоговые вкладки с параметрами расчета куба «por»

97

4.2. Расчет куба проницаемости

Скважинные данные о проницаемости представлены в колонке Rgi.Кпр. Особенностью проницаемости является ее резкая изменчивость. Такие данные трудно интерполировать. Гистограмму проницаемости по скважинным данным показывает синяя кривая на фиг. 95. Гистограмму проницаемости для куба после интерполяции скважинных данных (ее опции совпадают с опциями расчета пористости) показывает соседняя красная кривая. Можно видеть, что в результате интерполяции гистограмма проницаемости существенно искажается. Принять это искажение нельзя, даже делая скидку на неравномерное размещение скважин в кубе.

Фиг. 95. Гистограмма проницаемости по скважинным данным Rgi.Кпр (синяя кривая) и гистограмма того же параметра в кубе (красная кривая)

Попробуем выполнить интерполяцию логарифма проницаемости, который меняется не так резко. Логарифмирование кривых Rgi.Кпр производится при помощи калькулятора (фиг. 96).

98

Фиг. 96. Калькулятор кривых десятичного логарифма проницаемости Rgi.LogКпр

Гистограмма логарифма проницаемости по скважинным данным и гистограмма того же параметра после интерполяции в кубе показаны на фиг. 97. Результат интерполяции логарифма проницаемости можно считать приемлемым.

Фиг. 97. Гистограмма логарифма проницаемости по скважинным данным (синяя кривая) и гистограмма того же параметра в кубе (красная кривая)

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

99

показаны на фиг. 98. Результат интерполяции проницаемости через логарифм лучше предыдущего и признается удовлетворительным.

Фиг. 98. Гистограмма проницаемости по скважинным данным (синяя кривая) и гистограмма проницаемости, рассчитанной в кубе через логарифм (красная кривая)

Фиг. 99. Калькулятор, выполняющий обратное преобразование логарифма проницаемости в ячейках куба

100

4.3. Расчет куба нефтенасыщенности

Расчет куба нефтенасыщенности производится посредством интерполяции кривых нефтенасыщенности Rgi.Kн (фиг. 100). Главное отличие от интерполяции пористости состоит в том, что для Kн метод интерполяции устанавливается не «Согласно структуре», а «По горизонтали», поскольку нефтенасыщенность стратифицируется по вертикали. Еще одно отличие следует из того, что при интерполяции «По горизонтали» возможно появление зон коллектора, экранированных рельефом пласта и имеющих вследствие этого неопределенную насыщенность, как в неколлекторе. Для заполнения таких областей значениями Kн необходимо применять дополнительную корректирующую процедуру

«Tools\Interpretation\Заполнить по литологии».

Фиг. 100. Параметры расчета куба нефтенасыщенности

101

Фиг. 101. Сечение куба нефтенасыщенности целевого горизонта, до (вверху) и после (в центре) заполнения по литологии (внизу)

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

4.4. Построение карт подсчетных параметров

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

1.Карта эффективной нефтенасыщенной мощности коллектора;

2.Карта средневзвешенной пористости нефтенасыщенной части коллектора;

3.Карта средневзвешенной нефтенасыщенности нефтенасыщенной части коллектора.

Расчет карты эффективной нефтенасыщенной мощности коллектора производится по кубу NTG. В меню, ассоциированном с кубом NTG, выберем пункт «Средства\Calculate Map\Эффективная мощность». Откроется диалоговое окно расчета карты эффективной мощности (фиг. 102).

102

Фиг. 102. Расчет карты эффективной нефтенасыщенной мощности коллектора

Имя карты-результата «HefOil». Поскольку мощности считаются по кубу NTG, это гарантирует совпадение карты (при отсутствии сглаживания) с эффективными мощностями по скважинам. Поля «Top» и «Bottom» указывают, что расчет ведется между поверхностями «1» и «2». Поле «Bound 1» указывает, что расчет ведется только по части коллектора выше ВНК. Значение других опций диалогового окна описаны в «Руководстве пользователя», том 5, раздел 6. Например, если помимо ВНК мы имеем еще и ГНК, то поверхность ГНК помещается в поле «Bound 2», и в соседнем списке выбирается значение «Down».

Полученная карта эффективной нефтенасыщенной мощности показана на фиг. 103.

103

Фиг. 103. Карты эффективной нефтенасыщенной мощности коллектора

Расчет карты средневзвешенной (по вертикали) пористости для нефтенасыщенной части коллектора производится по кубу «por». В меню, ассоциированном с кубом «por», выберем пункт «Средства\Calculate Map\Средневзвешенные» Откроется диалоговое окно расчета карты средневзвешенного параметра (фиг. 104).

Фиг. 104. Расчет карты средневзвешенной пористости

Имя карты-результата «poro_oil». Имя исходного куба «por». Значения куба NTG используются как сомножители, поскольку куб «por» завышает пористость в ячейках, через которые проходит граница коллектор\неколлектор.

Чем средневзвешенные (по вертикали) значения параметра отличаются от средних (по вертикали) значений параметра? Это легко выяснить. Мы знаем, что все процедуры, вызываемые через пункт меню «Средства», написаны на языке Th. Соответствующие фрагменты кода показаны на фиг.105. При расчете средневзвешенных значений параметра vp весовыми коэффициентами являются мощности ячеек dz и величины v из куба NTG (фиг.105, вверху). При расчете средних значений параметра vp различия ячеек по толщине и куб NTG не учитываются (фиг.105, внизу). Фиг.105 иллюстрирует открытость DV-Geo в части кодов расчетных процедур.

104

Фиг. 105. Фрагменты исходного кода процедур расчета карт средневзвешенных (вверху) и средних (внизу) значений.

Полученная карта средневзвешенной пористости нефтенасыщенной части коллектора показана на фиг. 106.

105

Фиг. 106. Карта средневзвешенной пористости нефтенасыщенной части коллектора

Расчет карты средневзвешенной нефтенасыщенности для нефтенасыщенной части коллектора отличается от предыдущего только тем, что вместо куба «por» используется куб «nnas». Вызовем меню куба «nnas» и в этом меню выберем пункт «Средства\Calculate Map\Средневзвешенные». Полученная карта средневзвешенной нефтенасыщенности

«nnas_oil» показана на фиг. 107.

Фиг. 107. Карта средневзвешенной нефтенасыщенности нефтенасыщенной части коллектора

106

4.5. Построение контура ВНК

Все процедуры подсчета запасов используют в качестве полигона внешний контур ВНК. Когда подсчет запасов производится по зонам на плане XY (чисто газовой, газонефтяной, чисто нефтяной, водонефтяной), дополнительно требуются внешний и внутренний контуры ГНК, а также внутренний контур ВНК. Названные контуры должны быть представлены как объекты типа «Area». В качестве примера построим внешний контур ВНК. Для этого выполним следующие действия:

1.При помощи калькулятора рассчитаем разность между кровлей целевого горизонта (поверхностью «1») и уровнем ВНК (поверхностью «OWC»). Калькулятор показан на фиг. 108. Результатом будет поверхность «Difference», фиг. 109.

Фиг. 108. Расчет разности поверхностей «1» и «OWC»

2.Через меню поверхности «Difference» в «Списке объектов» создадим ее карту в изолиниях (объект «OWC»). В диалоговом окне расчета изолиний (открывается через меню карты в изолиниях «OWC») определим расчет единственной изолинии нулевого уровня. Это будет основная часть создаваемого контура ВНК. Бросим рассчитанную изолинию на план XY (фиг. 109). При расчете изолинии необходимо убедиться, что установленное значение «ноль» соответствует именно уровню, но не шагу изолиний на соседней закладке, иначе процедура закончится аварийно.

107

Фиг. 109. Расчет по карте «Difference» изолинии нулевого уровня

3.Преобразуем рассчитанную изолинию в объект типа «Area». Для этого захватим ее на плане XY сначала правой, а потом (удерживая правую кнопку) еще и левой кнопкой мыши. При этом рядом с курсором мыши появится рамка с именем «Полигон».

Перетащим изолинию в группу «Area\OWC» окна «Список объектов». В момент помещения изолинии в группу «Area\OWC» необходимо будет дать ей новое имя

(«Vnesh»).

4.Удалим с плана XY объект «изолиния» и бросим на план точно такой же по форме объект

«Area\OWC\Vnesh».

5.Контур ВНК (объект «Vnesh») должен быть замкнут. Замкнем его при помощи ручного редактирования вдоль линии разлома «Area\Raz\raz5_2_1». Для перехода в режим редактирования контура необходимо вызвать его диалоговое окно двойным нажатием левой кнопки мыши (с указанием на контур на плане XY), и отметить в этом окне элемент «Редактировать». Появится кнопочная панель «Редактировать полигон». Устанавливая на этой панели режимы «Редактировать» и «Замкнуть» (кнопки на противоположных краях верхнего ряда), контур «Vnesh» замыкается по линии разлома (фиг. 110). Построение контура ВНК завершается сохранением модели.

108

Фиг. 110. Замыкание контура ВНК по линии разлома

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

109

5. Подсчет запасов

5.1. Подсчет запасов по картам эффективных и средневзвешенных параметров

Для подсчета запасов по картам эффективных параметров необходимо через меню главного окна DV-Geo «Данные\Дополнительные средства» выбрать модуль «Инструменты\Подсчет запасов» (фиг. 111).

Фиг. 111. Вызов модуля «Подсчет запасов»

Главное окно модуля «Подсчет запасов» показано на фиг. 112. В этом окне необходимо ввести имя задания (допустим, «Primer_1»), выделить галочкой имя TH-Программы «Нефть и растворенный газ» и нажать кнопку «Создать задание».

Фиг. 112. Создание задания в модуле «Подсчет запасов»

В левой части окна появится пункт «Primer_1». В меню пункта «Primer_1» выберем команду «Добавить пласт», после чего введем имя «P7» - имя пласта в модели «CGE». Соответствующий пункт появится в левой части окна (фиг. 113).

110

Фиг. 113. Добавление пласта в модуле «Подсчет запасов»

Вышестоящий пункт «Дополнительные зоны» позволяет задать контуры лицензионного участка, 2-х километровой зоны и водоохранных зон. Мы оставляем данные поля пустыми.

Если в левой части окна выделить имя пласта «P7», то в правой части окна отобразится диалоговая форма, в которой необходимо указать имена трех карт – карты эффективной нефтенасыщенной мощности, карты средневзвешенной нефтенасыщенной пористости и карты средневзвешенной нефтенасыщенности (по нефтенасыщенной части коллектора, фиг. 114). Имена карт можно перенести мышью из окна «Список объектов».

Фиг. 114. Исходные карты в модуле «Подсчет запасов»

Если в левой части окна выделить имя зоны или категории (фиг. 114), то в правой части окна отобразится поле, в которое можно перенести один или несколько соответствующих контуров типа «Area», тоже мышью из списка объектов. Мы переносим сюда контуры зон ВНЗ (это внешний контур ВНК), ЧНЗ (это два внутренних контура ВНК), категории В (контур «b») и категории С1 (контуры «c1_1» и «c1_2»). Контуры категорий были построены геологом.

После заполнения всех названных полей следует вновь выделить в левой части пункт «Primer_1» и в соответствующей правой части окна нажать кнопку «Выполнить» (фиг. 113). Результатом работы процедуры будет таблица MS Excel с отображением промежуточных и окончательных результатов подсчета запасов, и с разбиением по категориям, и суммарных

(фиг. 115).

111

Фиг. 115. Результат подсчета запасов по картам эффективных и средневзвешенных параметров

5.2. Подсчет запасов с использованием кубов параметров

Для подсчета запасов непосредственно по кубам параметров необходимо через меню главного окна DV-Geo «Данные\Дополнительные средства» выбрать модуль «Инструменты\Оценка модели». Окно модуля показано на фиг. 116.

112

Фиг. 116. Окно модуля «Оценка модели»

Через меню корневого пункта «Оценка модели» в левой части окна вводится имя задания «Primer_2». При выделении имени задания в правой части окна модуля отображается диалоговое окно для ввода данных. Фиг. 116 показывает, как оно должно быть заполнено. С клавиатуры здесь вводится только название месторождения, объекта и залежи, остальные поля окна заполняются перенесением соответствующих элементов из «Списка объектов». В качестве полигона мы должны поместить рассчитанный ранее контур ВНК. В поле «Список скважин» можно поместить таблицу типа «Objects» (фиг. 117). Скважины (если их список представлен) используются для сравнения данных Rgi с параметрами кубов.

Фиг. 117. Таблица с именами скважин, помещаемая в поле «Список скважин»

113

Поясним, как можно быстро создать список скважин, отображаемых в окне «План горизонтального слайса». Для этого необходимо выполнить следующее. Первое - в меню плана XY (фиг. 29) выбрать пункт «В набор…», после чего откроется диалоговое окно «Набор скважин». Второе – захватить правой кнопкой мыши план XY и перенести его в центральное поле окна «Набор скважин». После этого в нем появятся номера скважин (фиг. 118). Третье – нажать кнопку «Выполнить». Процедура попросит дать имя для таблицы скважин, после чего набор данных под указанным именем появится в группе «Objects» окна «Список объектов».

Фиг. 118. Создание таблицы с именами скважин

Вернемся к заполнению окна модуля «Оценка модели». Выделяя в левой части окна пункты категорий B и C1, перенесем в списки соответствующих полигонов контуры «Area\kat\b» для B и «Area\kat\c1_1», «Area\kat\c1_2» для C1. Теперь все исходные данные в окне модуля определены. Нажатие кнопки «Выполнить» запускает расчетную процедуру, результатом которой является таблица MS Excel (фиг. 119).

114

Фиг. 119. Результат выполнения процедуры «Оценка модели»

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

Итак, мы выполнили подсчет запасов в нашей модели по двум методикам – по картам эффективных параметров и по кубам свойств. Что можно сказать о результатах, они совпали? Сравнительная таблица результатов показана на фиг. 120. Можно видеть, что нефтенасыщенные площади и объемы достаточно хорошо совпадают. Видимое различие по запасам нефти (когда во второй колонке берется число вне скобок) следует из того, что при расчете по кубам свойств не учитывается плотность нефти (0.949) и коэффициент приведения

115

к поверхностным условиям (0.972). Если учесть два названных коэффициента, то и по запасам совпадение будет хорошим (во второй колонке число в скобках).

Фиг. 120. Сравнение результатов подсчета запасов по двум методикам

Некоторые дополнительные подробности, касающиеся подсчета запасов, приведены в «Руководстве пользователя», том 5, раздел 8.

116

6. Подготовка данных для гидродинамического моделирования

6.1. Построение и редактирование сетки для гидродинамического моделирования

Стратиграфическая сетка геологической модели может включать миллионы ячеек. Это число слишком велико для расчетных гидродинамических систем (так называемых «симуляторов»). По этой причине последним этапом геологического моделирования обычно является ремасштабирование, или укрупнение сетки. Суть названного этапа состоит в том, чтобы при сокращении числа ячеек в десятки и сотни раз мы не потеряли существенные особенности модели.

Фиг. 121. Отображение ячеек стратиграфической сетки в режиме «Показать все ячейки»

Порядок создания гидродинамической сетки следующий. В меню группы «Lattice» в окне «Список объектов» выбираем пункт «Добавить сетку». Появляется диалоговое окно «Новая сетка», показанное на фиг. 122. В этом окне, помимо имени новой сетки, указываем имя куба, по которому она создается (из указанного куба берется только его геометрия), и начальное число узлов по каждому измерению. Необходимо также задать опции «Отрегулировать по кубу» и «Установить по поверхностям».

Выбор опции «Отрегулировать по кубу» ведет к тому, что границами сетки в пространстве 3D будут внешние границы куба. В противном случае сетка будет иметь контур куба только на плане XY, а сверху и снизу ее границами будут горизонтальные плоскости на уровне минимума и максимума куба по Z. Опция «Установить по поверхностям» в дополнение к предыдущей определяет, что разбиение на слои ячеек вдоль координаты Z производится лишь поверхностями структурного каркаса куба. В этом случае заданное выше количество узлов сетки по вертикали Nz игнорируется.

Созданную сетку можно визуализировать, перенеся ее из «Списка объектов» в окно «Model XYZ» или на план регулярного слайса. Для работы с гидродинамической сеткой используется кнопочная панель «Инструментарий сетки», которая открывается нажатием одноименной кнопки на панели главного окна DV-Geo (фиг. 123). Первое, что можно сделать при помощи

117

этой панели – это изменить процент отображаемых узлов сетки, так как по умолчанию она отображается только в режиме «Граница» (внешний контур).

Фиг. 122. Диалоговое окно «Новая сетка»

Фиг. 123. Инструментарий для работы с гидродинамической сеткой

118

Обратите внимание на следующую разницу. Хотя гидродинамическую сетку можно отобразить в окне обыкновенного регулярного слайса (как на фиг. 123), для ее редактирования необходимо открыть особое окно. Чтобы его открыть, надо выполнить следующие действия: первое - визуализировать сетку в окне «Model XYZ»; второе - создать нажатием кнопки на панели «Инструментарий сетки» слайс в кубе по линии узлов; третье - открыть нажатием соседней кнопки специализированное окно «слайс по линии узлов», в котором только и возможно названное редактирование.

Перемещать слайс по линии узлов можно мышью, захватывая его в кубе 3D. При этом должна быть нажата кнопка, разрешающая перемещение (справа от списка «Процент сетки»). Вторая возможность двигать слайс – при помощи кнопок пошагового перемещения на панели «Инструментарий сетки».

Покажем, как выполняется редактирование сетки. Добавление линий узлов производится следующим образом. Первое - на специализированной панели устанавливается режим «Выделить линию узлов» (фиг. 123). Второе - в окне редактирования сетки нажатием левой кнопки мыши выделяются две соседние однонаправленные линии узлов, между которыми мы хотим добавить одну или несколько новых линий. Третье - нажатием кнопки «Добавить линии узлов» открывается диалог, в котором мы указываем, на сколько частей мы хотим разбить интервал между выделенными линиями. После нажатия кнопки «OK» разбиение выполняется. Удаление выделенных линий производится нажатием кнопки «Удалить линии узлов». Есть режим смещения линий узлов, смещения отдельных узлов и другие. Основные режимы редактирования сетки можно освоить опытным путем, ориентируясь по «всплывающим» описаниям кнопок панели. Подробное описание работы с гидродинамической сеткой приведено в «Руководстве пользователя», том 2, раздел 3-8, и том 5, раздел 9.

Для того, чтобы усреднение свойств в ячейках гидродинамической сетки было минимальным, перед ее построением (и в первую очередь – перед построением горизонтальных линий узлов) строится геолого-статистический разрез (ГСР) по кубу «LitoPr». Порядок расчета названного разреза показан на фиг. 124. Кривая, отображающая разрез, помещается в группу «Stat», дочернюю для куба «LitoPr».

Фиг. 124. Расчет кривой геолого-статистического разреза

119

Кривая ГСР является объектом типа «Area». Для ее отображения необходимо открыть окно типа «Кросс-плот», что делается из меню главного окна DV-Geo (команда «Пространство\Кросс-плот\2D»). На фиг. 125 показана кривая ГСР куба «LitoPr» и ее использование при задании толщин ячеек гидродинамической сетки.

Фиг. 125. Задание толщин ячеек на основании геолого-статистического разреза

6.2. Пересчет кубов параметров на гидродинамическую сетку

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

«Upscale» (фиг. 126).

Фиг. 126. Вызов процедуры «Upscale» из меню гидродинамической сетки

Окно процедуры «Upscale» показано на фиг. 127. Поясним работу в этом окне.

120

Фиг. 127. Окно процедуры «Upscale» с установленными параметрами

Первое, что мы устанавливаем – мы выбираем из списка «Параметр пересчета» геологический куб, который мы собираемся ремасштабировать. Мы будем выполнять пересчет на гидродинамическую сетку куба пористости «Pack.por». Напомним, что саму сетку мы определили тем, что вызвали процедуру «Upscale» через ее меню.

Встроке над списком параметров пересчета задается имя результирующего куба на укрупненной сетке – мы задаем имя «por». В строке под тем же списком задаются числа дробления гидродинамической ячейки по направлениям, Nx, Ny, Nz. Это дробление не меняет сетку. Оно служит для отбора из геологического куба Nx*Ny*Nz. значений, равномерно распределенных по объему гидродинамической ячейки, для их последующего усреднения.

Вполе параметра «Точность» указывается доля заполнения гидродинамической ячейки значениями из геологического куба, которая должна быть превышена, чтобы ячейка считалась заполненной.

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

121

После выбора кубов в нижней части окна следует нажать кнопку «Добавить куб». В результате в верхней части окна появляются имена двух кубов на гидродинамической сетке, которые предстоит рассчитать. Их расчет производится последовательно, запуск расчета производится нажатием кнопки «Ремасштабирование». Результат ремасштабирования параметра «пористость» показан на фиг. 128. Можно видеть, что куб NTG (показывающий долю коллектора в ячейке) частично компенсирует потерю детальности вследствие увеличения размера ячеек.

Фиг. 128. Результат ремасштабирования пористости. Показан исходный куб пористости, ремасштабированный куб пористости, и ремасштабированный куб NTG

Для расчета куба проницаемости в том же окне используется закладка «A-H Average». Здесь можно выбрать направление и выполнить усреднение в пределах ячейки гидродинамической сетки в порядке «арифметическое-гармоническое» или наоборот. Арифметическоегармоническое усреднение означает, что сначала вдоль указанного направления производится арифметическое усреднение параметра, после чего все полученные значения усредняются гармонически. Напомним, что средним гармоническим x1, x2,…, xn является величина, обратная к среднему арифметическому 1/x1, 1/x2,…, 1/xn. Среднее гармоническое всегда меньше или равно среднему арифметическому, причем равно только в случае равенства

x1, x2, …, xn. Арифметическое-гармоническое среднее будет больше при выборе того направления, в котором параметр меняется сильнее.

После расчета гидродинамических кубов производится их экспорт в форматах систем «Eclipse» и «More». Процедура экспорта вызывается через меню кубов в списке объектов.

Дополнительные сведения о процедуре ремасштабирования смотрите в «Руководстве пользователя», том 5, раздел 9. Экспортом ремасштабированных кубов свойств геологическое моделирование в проекте «CGE» завершается.

122

Часть III. Геостатистическое моделирование в проекте «Primer»

Мы только что рассмотрели построение относительно простой детерминированной модели «CGE». В данном разделе мы будем говорить о построение модели на основе геостатистики. Напомним, геостатистикой называют группу методов интерполяции точечных данных, основанных на использовании вариограмм.

Начнем с простого вопроса: «Чем именно нас не устраивает детерминированная интерполяция?». Рассмотрим в этой связи сечение куба литологии, который мы получили для модели «CGE» (фиг.129). Этот куб точно соответствует имеющимся скважинным данным. Беда только в том, что его горизонтальная изменчивость в разбуренной части площади (левая половина сечения) сильно отличается от той же самой изменчивости в неразбуренной части площади (правая половина сечения). Понятно, что такой особенности в реальности быть не может.

Из картины на фиг. 129 можно сделать два вывода. Первый вывод (очевидный) состоит в том, что детерминированная интерполяция дает нам ложное представление о среде в правой части сечения, то есть там, где у нас недостаточная плотность данных. Второй вывод (менее очевидный) - мы не можем доверять всей полученной картине, пока мы не определим, какая плотность данных является достаточной.

К каким практическим негативным последствиям приводит детерминированное моделирование? К низкой точности прогноза? Нет, не только. Мы не можем говорить о низкой точности прогноза на правой части сечения на фиг. 129, когда картина принципиально неверна. Когда модель содержит столь очевидный изъян, наивно надеяться, что мы сумеем адаптировать ее к данным разработки.

Фиг.129. Сечение куба литологии модели «CGE»

123

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

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

Что надо будет сделать затем? Получив реализацию куба литологии, мы должны будем рассчитать для нее (в пределах коллектора) реализацию пористости и реализацию насыщенности. Так мы получим реализацию геологической модели. Желательно будет получить таких реализаций хотя бы десяток (а лучше несколько десятков). Для каждой из них мы должны будем подсчитать запасы. Поскольку реализации будут отличаться, в итоговых графах таблиц на фиг. 115, 119 мы получим некоторое распределение значений, объективно показывающее нам неопределенность подсчета запасов.

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

В настоящем учебном курсе мы пройдем весь этот путь от начала и до конца, то есть до оценки неопределенности объемов добычи. Когда мы действительно сумеем его пройти. Пока же мы сделаем только первый шаг – рассчитаем некоторое число реализаций кровли резервуара. Это будет сделано средствами DV-Geo на реальной модели, которую мы будем называть условным именем «Primer».

124

1. Построение реализаций кровли резервуара

1.1. Исходные данные

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

Фиг. 130. План размещения скважин и карта отражающего горизонта T по данным сейсморазведки 3D. Показана также линия сечения, представленного на фиг. 132

Наряду с картой по данным сейсморазведки 3D, имеются две карты той же поверхности, полученные в результате интерпретации сейсморазведки по сетям профилей 2D 1987 г. и 2000 г. Кроме того, имеется карта горизонта Т, полученная по данным гравиразведки (фиг. 131).

125

Фиг. 131. Карты горизонта Т, полученные по результатам интерпретации, слева направо, сейсморазведки 2D 1987 г., 2000 г. и гравиразведки

Вид представленных поверхностей горизонта Т в вертикальном сечении показан на фиг. 132. Линия сечения изображена на фиг. 130.

Фиг. 132. Поверхности T по результатам интерпретации: сейсморазведки 1987 г. (зеленая линия), сейсморазведки 2000 г. (голубая), сейсморазведки 2002 г. (красная), гравиразведки (серая линия). Черная линия - кровля доюрского основания по данным сейсморазведки 2002 г. Желтая линия - поверхность ВНК. На скважинах показаны отметки горизонта П2

Можно видеть, что данные сейсмических и гравиметрических работ разных лет покрывают различные участки площади и не совпадают друг с другом там, где их можно сравнить. Мы собираемся выполнить геостатистическую интерполяцию отметок П2 кровли резервуара на скважинах, используя синтезированную (на основе четырех имеющихся) поверхность горизонта Т в качестве тренда. Но на какой формальной основе мы можем объединить эти четыре поверхности? Чтобы ответить на этот вопрос, мы должны ранжировать имеющиеся поверхности по качеству. Для этого мы можем воспользоваться геостатистическим анализом.

126

1.2. Геостатистический анализ исходных данных

Суть геостатистического анализа будет в том, что мы исследуем, насколько хорошо каждая из четырех поверхностей (показанных на фиг. 130-131) коррелирует с отметками горизонта П2 на скважинах. Для этого мы будем придерживаться следующего плана действий:

1.Рассчитаем точки пересечения названных поверхностей со скважинами. Представим отметки поверхностей на скважинах как объекты типа «Correlation».

2.Рассчитаем разности отметок поверхностей и отметок кровли горизонта П2 на скважинах. Полученные разности (невязки) представим в виде таблиц типа «Objects».

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

4.Зная координаты пересечений поверхностей и скважин, рассчитаем вариограмму невязок по каждой таблице. Убедимся, что каждая вариограмма выходит на «плато» (что является еще одним признаком стационарности). Далее используем второй критерий качества поверхности

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

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

Приступим к выполнению нашего плана. Создание репера, фиксирующего точки пересечения поверхности со скважинами, производится через меню соответствующей поверхности в списке объектов (фиг. 133).

127

Фиг. 133. Расчет точек пересечения поверхности «T_DEPT_CGE_2002» со скважинами

Разности (невязки) отметок поверхности и отметок «P2_top» рассчитываются при помощи калькулятора (фиг. 134). Результат помещается в список объектов как таблица типа «Objects». Такая таблица рассчитывается для каждой из четырех исследуемых поверхностей.

Фиг. 134. Расчет невязок между ручным репером «P2_top» и отметками поверхности

«T_SEIS2002». Имя результата «DIF_P2_SEIS2002». Тип результата – «Data», то есть таблица типа «Objects»

Для расчета гистограмм и вариограмм невязок необходимо через меню главного окна DV-Geo «Данные\Дополнительные средства» выбрать модуль «Геостатистика\Геостатистический анализ». Окно модуля показано на фиг.135.

128

Фиг. 135. Окно модуля «Геостатистический анализ»

Первоначально правая часть окна пуста, а в левой части присутствует только один корневой пункт «Геостатистический анализ». Необходимо вызвать меню корневого пункта и выбрать в нем команду «Добавить объект 2D». Процедура предложит определить имя объекта 2D (мы указали «T1»), после чего в правой части окна открывается форма для введения его данных.

На фиг. 135 показано заполнение формы, необходимое для расчета всех нужных нам гистограмм и вариограмм. Мы будем менять в этой форме только имена таблиц невязок, последовательно перенося соответствующие таблицы из списка объектов. Поле «Список скважин» оставляется пустым, имена скважин берутся из таблиц невязок. В поле «Поверхность» обязательно должна быть указана некоторая поверхность, позволяющая определять координаты X,Y скважин. Названные координаты необходимы для расчета вариограмм. Мы поместили в это поле детерминированную поверхность кровли П2.

Через меню пункта «Объект 2D T1» последовательно добавляются гистограммы и вариограммы, пока только как дочерние пункты в древовидном списке. Для расчета той или иной гистограммы или вариограммы (по таблице невязок, указанной в форме объекта T1)

129

необходимо выделить соответствующий пункт курсором, заполнить ассоциированную форму и нажать на кнопку «Вычислить». Формы для расчета всех гистограмм и вариограмм показаны на фиг. 136 и 137. Для отображения рассчитанных гистограмм и вариограмм необходимо открыть соответствующие окна (одно для гистограмм, и одно для вариограмм) через меню корневого пункта «Геостатистический анализ». Отображаются те гистограммы (вариограммы), которые выделены галочками рядом с их именами.

Фиг. 136. Заполнение формы для расчета гистограммы

Фиг. 137. Заполнение формы для расчета вариограммы

130

Фиг.138. Гистограммы невязок относительно поверхности T по результатам интерпретации: сейсморазведки 1987 г. (зеленая линия), сейсморазведки 2000 г. (голубая линия), сейсморазведки 2002 г. (красная линия), гравиразведки (серая линия).

Гистограммы, полученные по четырем таблицам невязок, показаны на фиг. 138.

На что здесь надо обратить внимание? Первое - похожи ли эти гистограммы на распределение Гаусса. Второе – чему равно среднее значение невязок. Это среднее указывает, на сколько метров отметки горизонта П2 выше соответствующей поверхности. Третье – чему равна дисперсия невязок. Как было сказано выше, величина дисперсии показывает, насколько тесно поверхность связана с ручными отметками кровли П2. По критерию дисперсии невязок поверхности T распределились следующим образом – поверхность сейсморазведки 3D 2002 г. является лучшей, следующей идет поверхность сейсморазведки 2D 2000 г. и далее идет поверхность гравиразведки. Поверхность сейсморазведки 2D 1987 г. вызывает сомнения сразу по всем признакам – несимметричная форма распределения, слишком большое среднее отклонение и слишком большая дисперсия.

131

Фиг.139. Вариограммы невязок относительно поверхности T по результатам интерпретации: сейсморазведки 1987 г. (зеленая линия), сейсморазведки 2000 г. (голубая линия), сейсморазведки 2002 г. (красная линия), гравиразведки (серая линия). Показана также используемая далее модель одной из вариограмм.

Вариограммы, полученные по четырем таблицам невязок, показаны на фиг. 139.

Можно видеть, что три из четырех вариограмм выходят на «плато», и что пороги вариограмм невязок с хорошей точностью равны дисперсиям невязок. По критерию радиуса вариограмм невязок поверхности T в порядке убывания качества распределились в том же порядке, что и по критерию дисперсии. Лучшей является поверхность по данным сейсморазведки 3D 2002 г., затем поверхность по данным сейсморазведки 2D 2000 г., затем поверхность по данным гравиразведки. Поверхность по данным сейсморазведки 2D 1987 г. вызывает сомнения. Судя по вариограмме невязок, имеет место их нестационарность. Другими словами, поверхность 1987 г. и ручные отметки П2 систематически расходятся. Этим же объясняется и слишком большое среднее невязок, и их слишком большая дисперсия.

У лучшей поверхности по данным сейсморазведки 3D есть только один недостаток – она не покрывает всю интересующую нас область. То же относится и ко второй по качеству поверхности по данным сейсморазведки 2D. Поэтому имеющиеся карты поверхности Т были объединены следующим образом. В контуре полигона 3D бралась карта по данным сейсморазведки 3D. Вне названного контура бралась карта по данным сейсморазведки 2D 2000 г. Наконец, вне двух первых контуров бралась карта по данным гравиразведки. Для сглаживания стыков первый и второй контур были окружены переходными зонами шириной 2000 м. Процедура расчета комбинированной карты была написана на языке Th. Карта по данным сейсморазведки 2D (1987 г.) не использовалась.

132

1.3. Расчет реализаций кровли резервуара

Полученная комбинированная поверхность рассматривалась как детерминированный тренд при расчете реализаций кровли горизонта П2 по отметкам на скважинах. В ходе расчета реализаций использовалась сферическая модель вариограммы с порогом 25 единиц и радиусом 4000 м (показана фиолетовым цветом на фиг.139). Это некоторая натяжка по отношению к периферийным частям кровли резервуара, поскольку для них вариограмма невязок имеет меньший радиус и больший порог.

Фиг. 140. Окно модуля «Кригинг» и заполненная форма расчета реализаций кровли П2

Расчета реализаций кровли производился при помощи дополнительного модуля DV-Geo «Кригинг». Вызов модуля выполняется через пункт «Данные/Дополнительные средства» в меню главного окна системы. Окно модуля показано на фиг.140. Первоначально в левой части окна присутствует только один корневой пункт – «Кригинг». Через его меню мы должны

133

создать его дочерний объект «2D симуляции». Последнему мы дали имя «P2». Выделение вновь созданного объекта курсором ведет к отображению ассоциированной с ним диалоговой формы в правой части окна. В заполненном виде форма расчета реализаций кровли П2 показана на фиг. 140.

На что здесь необходимо обратить внимание? Основные три поля («Поверхность», «Поверхность внешнего дрейфа» и «Репер») заполняются переносом соответствующих объектов из окна «Список объектов». Мы предварительно создали на регионе «seis» новую поверхность «SIM_P2» (без реального расчета) и перенесли это имя в поле «Поверхность». Имя «T_SEIS3D2DGR» имеет комбинированная трендовая поверхность, созданная в предыдущем разделе. Имя «P2_top» имеет репер кровли резервуара на скважинах. Кроме этого устанавливается число реализаций, которое мы собираемся рассчитать (от 0 до 9, всего десять) и радиус окрестности кригинга (5000 м). Остальные параметры в полях формы задаются «по умолчанию». Тем не менее, надо проследить, чтобы они имели именно такие значения, как на иллюстрации.

Что еще необходимо для расчета реализаций? Осталась еще одна важная деталь - необходимо задать вариограмму. Вариограмма добавляется через меню пункта «Вариограммы», дочернего для пункта «P2 (2D Симуляция)» (фиг. 140). В названном меню необходимо выбрать команду «Новая вариограмма» и указать ее имя (мы указали «P2»). Форма для задания вариограммы и окно отображения вариограммы показаны на фиг. 141.

Фиг. 141. Форма для задания вариограммы и окно отображения вариограммы

Вариограмма для расчета реализаций задается в разделе формы «Модель вариограммы». Мы указали параметры модели, которую мы подобрали ранее на фиг. 139 (сферическая, радиус 4000, порог 25). Для расчета модели необходимо нажать в этом же разделе кнопку «Вычислить». Окно отображения вариограммы открывается нажатием кнопки «Показать».

134

Границы области отображения в окне вариограммы корректируются при помощи меню «Определить границы» и диалога «Границы координат».

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

Всего мы рассчитали 22 реализации кровли П2. Две из них показаны в окне «Профиль» на фиг. 142. Можно видеть, что на скважинах каждая реализация точно проходит через заданный репер. Между скважинами реализации следуют рельефу направляющей поверхности, рассчитанной по данным сейсморазведки 3D, 2D и гравиразведки. Отличия реализаций друг от друга соответствуют вариограмме невязок репера относительно тренда и показывают нам имеющуюся неопределенность.

Фиг. 142. Две реализации кровли резервуара в окне «Профиль». Показана также комбинированная поверхность, использованная в качестве тренда

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

135

Фиг. 143. Две реализации кровли резервуара вместе с уровнем ВНК (синяя поверхность) и кровлей ДЮК (доюрского комплекса, коричневая поверхность)

1.4. Оценка неопределенности объема нефтенасыщенных пород

Наличие 22-х реализаций кровли резервуара позволяет нам судить о неопределенности объема нефтенасыщенных пород. Это объем выше уровня ВНК и выше ДЮК, но ниже кровли П2. Расчет указанного объема выполняется для каждой реализации кровли П2 при помощи процедуры на языке Th. Все необходимые нам поверхности заданы на совпадающих сетках 50м*50м, что сильно упрощает вычисления. Мы производим перебор ячеек сетки, для каждой ячейки находим более высокую поверхность из пары (ВНК, кровля ДЮК), и определяем расстояние h (в метрах) от нее до П2. Если поверхность П2 выше, в искомый объем

136

добавляется 50*50*h кубических метров. Дополнительное ограничение состоит в том, что расчет производится только в контуре лицензионного участка (центральное поднятие на фиг. 143), поднятия на окраинах исключаются.

Полученное распределение значений объема нефтенасыщенных пород (миллионов кубических метров) по 22 реализациям показано на фиг. 144. Оно похоже на гауссовское. Среднее значение – 3862 миллионов м3. Среднеквадратичное отклонение σ – 96 миллионов м3. Как известно, величина, распределенная по Гауссу, с вероятностью 95% попадает в интервал «среднее плюс-минус два σ», то есть объем нефтенасыщенных пород с названной вероятностью равен 3862 ± 192 миллиона м3.

Фиг. 144. Распределение объема нефтенасыщенных пород по 22 реализациям кровли П2

Точность получилась очень высокой – максимальная ошибка (два значения σ) составляет всего 5% от прогнозируемой величины. Однако надо иметь в виду, что мы решили лишь первую и самую простую часть задачи – геометрическую. Конечная точность подсчета запасов будет существенно ниже, поскольку к геометрической неопределенности добавится неопределенность прогноза литологии, пористости и насыщенности.