Федеральное агенТство по образованию
ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«ТЮМЕНСКИЙ ГОСУДАРСТВЕННЫЙ НЕФТЕГАЗОВЫЙ УНИВЕРСИТЕТ»
Методические указания
к дипломным и курсовым работам и
к лабораторным работам по дисциплине
«Основы компьютерных технологий решения геологических задач»
для студентов очной и заочной формы обучения специальности 130304 «Геология нефти и газа»
Тюмень 2006
Утверждено редакционно-издательским советом
Тюменского государственного нефтегазового университета
Авторы: Белкина В.А. - к.ф.-м.н., доцент
Василевская М.А. - ассистент
Ó Тюменский государственный нефтегазовый университет
2006
ВВЕДЕНИЕ
Компьютерные технологии - неотъемлемая часть интерпретации геологических измерений. В процессе интерпретации на разных этапах необходимо построение геологических карт (структурных, общих, эффективных, эффективных нефтенасыщенных толщин, пористости, проницаемости и т.д.). В современных пакетах геологического моделирования, как правило, существует большое число алгоритмов построения карт. Обоснованный выбор алгоритма, учитывающего особенности отстраиваемой карты (например, её дифференциальные свойства или статические характеристики), позволяет заметно повысить точность моделей при заданном наборе исходных данных. Важность этого обстоятельства в геологии обусловлена тремя моментами:
1. ограниченностью эмпирических данных;
2. наличием ряда ошибок в измерениях (случайных, систематических и т.д.);
3. высокой стоимостью полевых измерений.
Из сказанного видно, что специалист, занимающийся интерпретацией геологических данных, должен разбираться в методах моделирования, причем не только на геологическом уровне, но и на математическом. Поэтому в курсах «Основы компьютерных технологий решения геологических задач» и «Системы построения и мониторинга эксплуатационных моделей нефтегазовых объектов» предусмотрено освоение методов моделирования различных геологических параметров. Изучение алгоритмов картирования ведется на примере программного пакета Surfer (версия 8.0). В указании излагаются основные особенности и приёмы работы с Surfer 8.0, а также рассмотрены алгоритмы построения карт эффективных нефтенасыщенных толщин на разных стадиях разведки.
1. ПРОГРАММНЫЙ ПАКЕТ SURFER
1.1. Основные возможности пакета Surfer
Surfer позволяет строить и визуализировать двумерные карты, которые математически описываются функцией вида z=f(x,y). Назовем три основные функции пакета:
а) построение цифровой модели поверхности;
б) вспомогательные операции с цифровыми моделями;
в) визуализация поверхности.
Постановка задачи построения геологической карты формулируется как переход от набора значений функции Z в произвольных (неупорядоченных) точках плоскости (n точек с координатами Xi,Yi, i=1,2,…, n) к значениям этой функции в узлах некоторой регулярной сетки. В более общей постановке эта задача сводится к вычислению значений функции в любой точке области картирования (а значит, и в узлах сетки) по исходному набору данных.
Переход к значениям в регулярной сетке необходим для решения двух задач:
1) построения непрерывной поверхности (например, в виде карты), по которой можно было бы узнать значения параметра в любой точке области построения карты;
2) представления этих данных в виде математической цифровой модели, с помощью которой, решают другие геологические задачи, например, прогноз изменения геолого-промысловых параметров в межскважинном пространстве или подсчет запасов.
Пакет Surfer предоставляет широкий набор математических методов интерполяции и аппроксимации, которые позволяют выбрать наиболее подходящий алгоритм для повышения точности модели. Достоинством пакета Surfer является возможность проводить оценку качества исходных данных и получаемых результатов.
Перечисленные возможности пакета Surfer 8.0 позволяют проводить выявление "подозрительных" данных, и некачественных замеров. Имеется возможность определять участки, где наблюдается недостаток полевых данных. Таким образом, речь идет уже о задаче выбора оптимальной конфигурации размещения сетки скважин, обеспечивающую необходимую точность геологической модели. В пакете Surfer решение подобных задач обеспечивается построением карт значений погрешностей интерполяции, а также расширенными возможностями моделирования с использованием вариограмм.
1.2. Подходы к построению карт
Существует два основных принципа построения карт: в первом случае используется система треугольников, где вершины треугольников находятся в точках замера, во втором - регулярная сетка, где точки размещаются в узлах регулярной сетки (рис. 1.1).
При способе треугольников точки соседних скважин соединяют на плане линиями таким образом, что образуется система треугольников (рис. 1.1 а). Затем на каждой линии по правилу линейной интерполяции находят точки со значениями абсолютных отметок, кратными выбранной величине сечения между изогипсами.
Линейная интерполяция предполагает, что наклон линии, соединяющий две скважины, на всем её протяжении постоянен. Расстояние любой изогипсы от одной из точек наблюдения на этой линии при линейной интерполяции можно найти по формуле
Lx=[(Hх- H1)/( H2- H1)] l1,2 , (1.1)
где Lx - расстояние от искомой изогипсы до скважины 1 на линии, соединяющей скважины 1 и 2; Hх – значение (абсолютная отметка) искомой изогипсы; H1 и H2 – значения картируемой поверхности соответственно в скв. 1 и 2; l1,2 - расстояние между скв. 1 и 2.
Чем больше точек наблюдения, тем меньше размер треугольников и тем точнее построенная карта будет отражать форму реальной картируемой поверхности.
Одним из достоинств метода триангуляции является то, что при достаточном количестве экспериментальных точек по обе стороны от линии разрыва, цифровая функция отобразит этот разрыв.
. | |||||||||||
. | |||||||||||
. | . | ||||||||||
. | |||||||||||
. | |||||||||||
. | . | ||||||||||
. | |||||||||||
а) б)
Рис. 1.1. Основные принципы построения карт:
а) триангуляция; б) сеточный метод
Grid – представление цифровой модели значениями параметра в узлах регулярной сетки.
Интерполяция это один из способов аппроксимации. Аппроксимация (от лат. approximo — приближаюсь), замена одних математических объектов другими, в том или ином смысле близкими к исходным. Она позволяет исследовать числовые характеристики и качественные свойства объекта, сводя задачу к изучению более простых или более удобных объектов (например, таких, характеристики которых легко вычисляются или свойства которых уже известны). Интерполяция - это восстановление функции в промежуточной точке по известным её значениям в соседних точках, эти точки Хi и Yi, i=0,1,2,…,n называют узлами интерполяции, а функцию – интерполирующей или интерполянтом. Вид функции определяет способ интерполяции. На практике в качестве интерполирующей функции часто используются алгебраические полиномы различного порядка, так как полиномы легко вычислять, дифференцировать и интегрировать. Эта интерполяция называется полиномиальной. В простейшем (одномерном) случае задача интерполяции состоит в том, что по заданным точкам (Хi, Yi, i=0,1,2,…, n), требуется найти функцию F(Х), которая проходит через эти точки, т.е. выполняются равенства F(Хi)= Yi, i=0,1,2,…, n.
Все методы интерполяции делятся на 2 группы: детерминированные и геостатистические. Детерминированные методы для интерполяции используют математические функции (зависимости). Геостатистические методы базируются на математических и на статистических функциях, которые могут быть использованы для построения поверхностей и оценки точности прогнозов.
Многие карты отличаются большой сложностью, а именно сильной изменчивостью, разрывностью функций и т.д. При построении таких карт в целом по всей области [X1,Xn]*[Y1,Yn] необходимо использовать довольно сложные функции с большим числом параметров. Например, при интерполяции алгебраическими полиномами необходимо рассматривать полиномы высоких степеней. Но как известно из теории интерполяции, такие полиномы дают большие вычислительные погрешности. С другой стороны, при интерполяции «в целом» по всей области, при построении интерполирующего полинома используются значения Zi, расположенные друг от друга на значительном расстоянии, из чего может следовать их статистическая независимость.