Технология построения карт

 

Технология построения карт

 

Построение всех представленных в отчёте карт было выполнено в пакете «Cube-Technology» по результатам слежения временных разрезов и параметрам, полученным при ручной и автоматизированной корреляции данных ГИС (технологии ПИК). Ручная корреляция данных ГИС выполнялась по исходным каротажным кривым (PS), а автоматизированная по исходным кривым (PS) и результатами их обработки выполненной СК «Петроальянс» (относительной амплитуды , пористости, глинистости, признака насыщения).

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

В начале была выполнена кинематическая и динамическая увязка временных разрезов партий 16\88 и 08\91 перекрывающих Южно Ягунское месторождение. Увязка произведена следующим образом:

·         Уравнивание амплитуд введённых разрезов путём потрасной нормировки в скользящем окне 500 мс.

·         Горизонтальный спектральный анализ разрезов, выбор общей полосы частот 15-60 Гц. для всех разрезов, полосовая фильтрация в этой полосе.

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

·         Фазовая коррекция временных разрезов по информации на пересечении профилей (переворот фаз, ротация фаз)

·         Прослеживание наиболее интенсивного отражающего горизонта кровли баженовской свиты на всех разрезах

·         Приведение всех разрезов к одной линии приведения путём минимизации невязок интенсивного отражения и ввода постоянных поправок минимизации во временные разрезы

·         Картирование прослеженного интенсивного отражения

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

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

Следующим этапом структурных построений была стратиграфическая привязка отражений к пластам исследуемой толщи БС10-БС11. Для этого рассчитаны глубинные разрезы для пяти профилей. Глубинные разрезы рассчитывались путём миграции Кирхгофа со скоростями, вычисленными по сейсмокаротажу разведочной скв. 158. Миграция выполнена в пределах зоны Френеля 0.005 с. К рассчитанным глубинным разрезам были привязаны каротажные кривые PS всех ближайших разведочных и эксплуатационных скважин. Привязка проведена на основании разности измеренных и абсолютных отметок пласта БС10-1 по данным СК «Петроальянс». На основании результатов привязки на глубинных разрезах с точностью до фазы были стратифицированы кровли пластов БС10-1 и кровли родниковой пачки 1. По волновой картине хорошо видно, что отражение в районе кровле пласта БС10-1 на глубинных разрезах соответствует интенсивной красной фазе на времени 1.94 сек. на временных разрезах.

Выбранная фаза была уверенно прослежена на всех временных разрезах. При слежении не возникло никаких неопределённостей при переходе от одного профиля к другому по любой композиции профилей, что подтверждает достоверность прослеженного горизонта. Затем, для кровли пласта БС10-1 была рассчитана карта t0 по следующей методике:

·         Скользящее сглаживание линий t0 на базе 200 м.

·         Минимизация невязок на пересечении профилей путём нахождения оптимальных постоянных составляющих линий t0 на профилях

·         Разброс остаточных невязок на пересечениях профилей на базе 2400 м.

·         Разрядка пикетов на профилях до шага 100 м.

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

·         Площадное сглаживание регулярной сетки на базе 200 на 200 м. с подгонкой значений в регулярной сетки к значениям t0 в пикетах профилей[1].

 

По этой же методике были построены карты t0 подошвы родниковой клиноформы и кровли баженовской свиты. Все построенные карты основаны только на сейсмических данных. Для построения аналогичных структурных карт необходимо использование сейсмокаротажа или отметок глубин в скважинах. Абсолютные отметки кровли пласта БС 10-1 приведённые «ПетроАльянсом» по ряду скважин определены с ошибками. Это хорошо видно на карте абсолютных отметок, где отдельные скважины создают не правдоподобные аномалии (рис 6.3).

 

Ошибки в абсолютных отметках подтверждаются также низким коэффициентом корреляции абсолютных отметок и значений t0 в скважинах вычисленных путём площадной линейной интерполяции t0 снятых на профилях в точки скважин. Он равен всего лишь 0.42. В то же время коэффициент корреляции отметок глубин найденных на основе этих t0 по данным сейсмокаратажа составляет 0.71. Основываясь на том, что глубины по сейсмокаротажу ведут себя аналогично t0, мы нашли способ отбраковки эксплуатационных скважин с аномальными абсолютными отметками. Он заключается в следующем: строится структурная карта* по t0 с учётом отметок глубин в скважинах по данным сейсмокаротажа. С полученной таким образом структурной карты снимаются значения глубин во всех эксплуатационных скважинах. Среднеквадратичная ошибка глубин по отношению к абсолютным отметкам в этих скважинах, составила 15 м. После чего были выбраны лишь те скважины, разность глубин в которых между абсолютными отметками и структурной картой не превосходит 5 м. Теперь строится окончательная структурная карта с учётом лишь таким образом выбранных скважин. Коэффициент корреляции абсолютных отметок глубин выбранных скважин с окончательной структурной картой составил теперь уже 0.67. Заметим, что другой принцип отбраковки скважин по большому углу наклона не привёл к повышению коэффициента корреляции.

На рис 6.4 приводится окончательная структурная карта кровли пласта БС10-1 построенная по сейсмическим данным с учётом абсолютных отметок в скважинах оставшихся после отбраковки.

В этом построении использованы также и разведочные скважины коэффициент корреляции абсолютных отметок, которых с t0 составляет 0.6. Построенная структурная карта проходит через отметки глубин в выбранных эксплуатационных и разведочных скважинах с точностью до 1 м. На рис. 6.5 приводится структурная карта кровли баженовской свиты. Она построена по рассмотренной в сноске методике с использованием отметок глубин определённых по сейсмокаротажу. Эта структурная карта имеет коэффициент корреляции с отметками глубин по сейсмокаротажу равный 0.96, что подтверждает её высокую достоверность. Структурная карта по подошвы родниковой клиноформы построена аналогичным способом приведена на рис 6.6 По структурной карте подошвы родниковой клиноформы рассчитана карта модуля градиента (максимального наклона) этой поверхности (рис 6.7)

 

 

 

На рис 6.8-6.10 приводится карты разности построенных структурных карт.

При построении карт t0 и структурных карт мы уделили особое внимание точности построений в условиях имеющейся у нас сейсмической и скважинной информации по Юго-Западному поднятию Южно Ягунскому месторождению и прилегающей территории. Однако даже при расчётах, подтверждаемых статистически, мы не можем с полной уверенностью идентифицировать выявленные на структурных картах аномалии. Это связано с тем, что по данным, изложенным в отчёте [1] верхняя часть разреза территории Южно Ягунского месторождения осложнена зонами вечной мерзлоты: «Южно-Ягунское месторождение приурочено к южной геокриологической зоне, характеризующейся двухслойным строением (1-й слой современной мерзлоты и 2-ой слой древней мерзлоты). Многолетнемерзлые породы первого слоя, как правило, приурочены к массивам торфяников. Кровля верхнего слоя мерзлоты (слоя современной мерзлоты) имеет четкую зависимость залегания от поверхностных условий и литологического состава грунтов. Подошва слоя современной мерзлоты в большинстве случаев лежит на глубинах 15 – 30 м с экстремальными значениями 50 – 80 м. Слой древней мерзлоты (С.Д.М.) имеет прерывистое распространение. Его отсутствие фиксируется под наиболее крупными речными долинами». Можно лишь надеяться, что карты разностей структурных карт не подвержены влиянию зон вечной мерзлоты, так как они получены путём вычитания глубин одного горизонта из другого, в результате чего влияния верхней части разреза, осложненное пространственными низкочастотными аномалиями, вызванными вечной мерзлотой и зонами её растепления, устраняется.

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

[1] Суть итерационного алгоритма площадного сглаживание регулярной сетки с учётом исходных значений в нерегулярной сетке состоит в следующем: пусть  множество значений функции в некоторой совокупности точек. С помощью линейной интерполяции можно вычислить множество значений в другой совокупности точек . Обратной линейной интерполяцией  нельзя получить исходные значения :

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

Начальное значение .

Иногда бывает достаточно одной итерации:

В этом случае уточненная интерполяция вычисляется по формуле:

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

 

* Методика построения структурной карты с учётом отметок глубин в скважинах

состоит в следующем:

 

Ø         расчет значений t0 в точках скважин по регулярной сетке t0

Ø         расчет единого линейного вертикального годографа Vср (t0) путём линейной аппроксимации Vср =A+B•t0 значений Vср =2H/t0 и t0 определённых в скважинах (H – отметки глубин в скважинах)

Ø         расчёт постоянной составляющей A в каждой скважине по формуле A = Vср -B•t0,

Ø         расчет регулярной сетки A (так же как для Vср )

Ø         расчет регулярной сетки Vср =A+B•t0 по регулярном сеткам A и t0

Ø         расчет глубин структурной карты по регулярным сеткам t0 и Vср