На диске D в папке Student находится 25 папок с именами «Вариант 1», «Вариант 2»,…, «Вариант 25», в каждой из них файл Пласт_Б10.dat, который нужно скопировать в свою папку. Этот файл содержит следующие параметры (столбцы):
1. Координата Х, м,
2. Координата Y, м,
3. Альтитуда А, м,
Рис. 3.8. Пример оформления геологической карты
4. Удлинение L, м,
5. Глубина кровли пласта Н, м,
6. Толщина общая h об, м,
7. Толщина эффективная h эф, м,
8. Толщина эффективная нефтенасыщенная h нэф, м,
9. Абсолютная отметка ВНК, м,
10. Открытая пористость m o, д. е.,
11. Наименование скважин.
Задание 1. Тема: «Сравнение методов построения цифровых геологических моделей». Построить структурную карту по кровле пласта Б10 восемью методами: крайгинга, радиальных базасных функций, обратных расстояний, минимальной кривизны, полиномиальной регрессии, триангуляцией с линейной интерполяцией, ближайшего соседа.
Замечание. Построение структурных карт ведётся по абсолютным отметкам (АО). В файле с исходными данными содержатся значения глубин, альтитуд и удлинений. Для вычисления абсолютных отметок воспользоваться выражением:
АО = -(Н-А-L). (4.1)
Структурную карту по кровле, построенную методом Крайгинга, оформить: нанести скважины, написать заголовок, вынести шкалу цветов, над шкалой цветов написать «Условные обозначения», под этой строкой указать единицу измерения. Вывести одновременно на экран все карты, сделать визуальный анализ карт и написать вывод, исходя из предположения, что залежь структурного типа вытянутая в меридиональном направлении. В выводе отметить, какие методы позволили построить структурную карту по кровле пласта, соответствующую геологическим представлениям.
Задание 2. Тема: «Изучение геометрии залежи». Построить карту общих толщин (hоб) и карту по подошве пласта Б10. последнюю рекомендуется получить как сумму цифровых моделей по кровле пласта и hоб. Вычислить среднее значение ВНК. В изучаемой модели залежи поверхность ВНК предполагается горизонтальной.
Получить в цифровом виде внутренний (по структурной карте подошвы) и внешний (по структурной карте кровли) контуры нефтеносности и нанести их в соответствии с принятыми в структурной геологии обозначениями на карту кровли. Полученные файлы с контурами нефтеносности назвать: vnutr.bln и vneshn.bln.
Вывести на экран все карты, характеризующие геометрию залежи. Они должны быть оформлены в соответствии c требованиями задания 1. Написать вывод. В выводе указать:
· перечень карт, описывающий геометрию залежи;
· описать тип залежи, изменчивость общей толщины на основе построенных карт.
Задание 3. Тема: «Построение карты эффективных нефтенасыщенных толщин h нэф с использованием априорной информации (значений h нэф, равных нулю на внешнем контуре нефтеносности)». Построить две карты h нэф:
1. только по замеренным значениям - h 1нэф;
2. по замеренным значениям и значениям нулей в ряде точек внешнего контура нефтеносности – h 2нэф.
Для построения карты h 2нэф следует в файл Пласт_Б10.dat из файла vneshn.bln скопировать столбцы координат Х и Y, начиная со второй строки. В столбец h 2нэф в добавленные строки занести нули. Для выполнения этих действий необходимо открыть два окна Worksheet.
Вывести на экран одновременно обе оформленные карты h 1нэф и h 2нэф. Написать вывод, в котором отметить:
· чем принципиально отличаются исходные данные карт;
· какая из карт отвечает геологическим закономерностям об изменении эффективных нефтенасыщенных толщин и почему.
Задание 4. Тема: «Построение карты эффективных нефтенасыщенных толщин с учетом геометрии залежи». Этот алгоритм основан на закономерности:
h нэф = h ноб * k песч, (4.2)
где h ноб – общая толщина нефтенасыщенной части залежи, м и k песч – коэффициент песчанистости, д. е.
Отметим, что карта h ноб вначале строится отдельно для чисто нефтяной зоны (ЧНЗ) и для водонефтяной зоны (ВНЗ). Обозначим их соответственно: h обчнз и h обвнз. Раздельное построение карт связано с тем, что нижней границей залежи в ЧНЗ и ВНЗ являются разные поверхности (Рис. 4.1): в ЧНЗ – структурная карта по подошве, а в ВНЗ – поверхность водонефтяного контакта.
Рис. 4.1. Геологический разрез залежи структурного типа
Схема возможного расположения чисто нефтяной и водо-нефтяной зон показана на рис. 4.2.
а б
Рис 4.2. Схема: а) ЧНЗ и б) ВНЗ
Рис. 4.3. Расположение зон ЧНЗ и ВНЗ в плане.
Окончательно карта общих толщин нефтенасыщенной части залежи получается по выражению:
h ноб = h обчнз + h обвнз. (4.3)
В данной работе k песч находится по формуле:
, i=1,2,…,n. (4.4)
Значение k песч для формулы 4.2 находится как средне-арифметическое выборки, полученной по формуле 4.4.
Итак, схема алгоритма такова:
1. Построить карту общих толщин h 1нэф как разность структурных карт подошвы и кровли.
2. Т. к. карта h 1нэф для рассматриваемой задачи актуальна только в ЧНЗ (рис. 4.2а.), надо удалить значения из цифровой сеточной модели за пределами внутреннего контура нефтеносности. Достичь этого можно командой Grid/Blank; в качестве файла типа *.bln следует использовать файл vnutr.bln.
3. Полученная в п. 2 карта в отбланкованной части содержит код отсутствия значений. Недостатком пакета Surfer является невозможность использования такого рода карт в арифметических преобразованиях. Поэтому необходимо код отсутствия значений заменить на нули. Это можно сделать, например, арифметическим выражением С=МАХ(А, 0) в команде Grid/Math. Под А здесь понимается карта, полученная в п. 2 а С – имя выходного файла. Полученную карту обозначим h обчнз.
4. Для получения карты h обвнз вначале вычислить карту h 2нэф как разность ВНК (в рассматриваемой модели залежи это горизонтальная плоскость, значит, он задаётся константой) и структурной карты по кровле.
5. Так как значения h 2об имеют смысл только в ВНЗ, необходимо удалить значения из сеточной модели за внешним контуром нефтеносности и внутри ЧНЗ (рис. 4.2б). Это достигается командой Grid/Blank, которую необходимо выполнить дважды. При использовании файла vneshn.bln отбланковать область вне контура, а при задании файла vnutr.bln – внутри контура.
6. В полученной в п. 6 карты заменить коды отсутствия значений на нули (см. п.3) - h обвнз.
7. Построить карту общих толщин нефтенасыщенной части пласта по формуле 4.3.
8. Вычислить значения k iпесч, i=1,2,…,n по формуле (4.4). Выполняется это в окне Worksheet в файле типа *.dat. затем найти средне-арифметическое значение k песч командой Data/Statistics.
9. Построить карту h нэф по формуле (4.1).
10. Вывести на экран одновременно две карты h нэф из заданий 3 и 4 и написать вывод. В выводе:
· Описать алгоритм построения карты h нэф по геометрии залежи;
· Сравнить карты h нэф построенные разными алгоритмами. Найти различия между ними и объяснить причины этих различий.
Задание 5. Тема: «Построение карты линейных запасов (Q лин) нефти».
Под картой Q лин понимается карта, полученная по формуле:
Q лин = h нэф *m o *k нн, (4.5)
где m o – открытая пористость;
k нн – коэффициент нефтенасыщения.
Отметим, что m o и k нн может быть как карты, так и постоянные значения. В изучаемой модели залежи имеется ряд замеров m o, поэтому целесообразно построить карту открытой пористости. Коэффициент k нн представлен только одним значением k нн =0,45.
В этой работе требуется:
1. Построить карту Q лин по формуле (4.5). Нанести на неё скважины, внутренний и внешний контуры нефтеносности, заголовок и все необходимые условные обозначения;
2. Вычислить запасы нефти (Q бал) и площадь залежи S на основе карты линейных запасов;
3. Написать вывод, в котором указать Q бал и S, написать их единицы измерения и обязательно отметить, для каких условий оценены запасы. По карте Q лин выделить участки с наибольшей и наименьшей плотностью запасов.